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(D ■ ABSTRACT 

We describe the calibration and data processing methods used to generate full-sky maps of 
the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy 
Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based 
, largely on analyses of the flight data supplemented, where necessary, with results from ground 

I tests. The data are calibrated in flight using the dipole modulation of the CMB due to the ob- 

(NJ ■ servatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative 

, algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pix- 

' elized sky map temperatures. The noise properties are determined by analyzing the time-ordered 

■ data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter 

I to the time-ordered data to remove a low level of 1// noise. We infer and correct for a small 

pi i' (^1%) transmission imbalance between the two sky inputs to each differential radiometer, and 

Q . we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. 

| H I No other systematic error corrections are applied to the data. Calibration and baseline artifacts, 

^ ' including the response to environmental perturbations, are negligible. Systematic uncertainties 

are comparable to statistical uncertainties in the characterization of the beam response. Both 
are accounted for in the covariance matrix of the window function and are propagated to un- 
' certainties in the final power spectrum. We characterize the combined upper limits to residual 

^ , systematic uncertainties through the pixel covariance matrix. 
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Subject headings: cosmic microwave background, cosmology: observations, space vehicles: in- 
struments, instrumentation: detectors 

1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) has produced full-sky maps of the cosmic mi- 
crowave background (CMB) of unprecedented precision and accuracy. On angular scales larger than ~ 0?5, 
the dominant uncertainty is not the instrument noise but rather the "cosmic variance" inherent when ana- 
lyzing a single realization (the observable universe) of a random process. In the cosmic variance limit, no 
further improvement can be made by reducing instrument noise, placing even greater importance on the 
minimization of non-random instrumental effects in the data. 

The WMAP design emphasizes control of systematic errors (Bennett et al. 2003a). The observatory 
was designed with a detailed systematic error budget in place, and a mature data analysis pipeline was 
written early to help inform many of the design decisions. Differential radiometers compare the temperature 
from widely-separated regions of the sky through back-to-back matched optics. Common-mode signals thus 
cancel before affecting the sky maps. The radiometer feed horns only illuminate a fraction of the primary 
mirrors, reducing the sidelobe response in the beam patterns. The instrument was designed to have minimal 
response to electrical or thermal perturbations and operates in an exceptionally stable environment at the 
second Earth-Sun Lagrange point. The observatory's compound spin and precession allow rapid inter- 
comparison of different positions on the sky, greatly reducing the coupling of systematic error signals into 
the sky maps and effectively symmetrizing the beam response. WMAP data are calibrated in flight using 
the dipole modulation of the CMB from the observatory's orbital motion around the Sun as a full-beam 
calibration source. We measure the beam pattern in flight using observations of the planet Jupiter. 

We characterize or limit systematic errors in the WMAP flrst-year data using flight data supplemented 
where necessary with results from ground tests. Systematic errors may be classified into several broad 
categories including the following: 

Calibration Errors. — The time-ordered data is simultaneously fit for the calibration and sky map. An 

iterative algorithm updates the calibration solution based on the previous iteration of the sky map solution, 
and vice- versa. The most important source of error is confusion between the dipole signal and higher-order 
sky signal, especially bright Galactic foreground emission in the low frequency WMAP bands. See §2.2.1. 

Map-making errors. — These are due to poor convergence in the sky map solution, or to errors in the 
determination of the spacecraft pointing. See §2.1.1 and §3.3. 

Beam Errors. — Instrument noise, background subtraction, and pointing errors can limit the in-flight 
measurement of the beam response from Jupiter. Although the beams are not symmetric, the observatory's 
compound spin and precession effectively symmetrize the beam response. Uncertainties in both the beam 
solid angle and the window functions must be characterized. See Page et al. (2003a) for a complete discussion 
of beam mapping. We summarize and incorporate their results in §3.3.3. 

Sidelobe Response. — Sidelobe pickup of bright sources (e.g. the Galactic plane) introduces an additive 
signal dependent on the orientation of the beams on the sky. Barnes et al. (2003) discuss the sidelobe 
response of each radiometer and estimate the effect on the first-year sky maps. 
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Baseline Errors. — Thermal or electrical perturbations can produce signals dominated by an additive 
term in the time domain. Slow drifts are removed as part of the calibration procedure, but signals near the 
spin period can couple to the sky maps with some efficiency. See §3.4. 

Striping. — Correlations in the time-ordered data from sources not fixed on the sky (e.g. 1// noise or 
post-detection filtering) introduce correlated noise in the sky maps. Application of a pre-whitening filter to 
the time-ordered data reduces this effect. See §2.3.2. 

We have constructed a detailed model of the instrument that successfully reproduces all major aspects 
of the instrument performance. Software simulations using this model validate the map-making algorithm 
and allow us to assess the effect in the sky maps of various signals in the time-ordered data. Based on this, we 
apply a pre-whitening filter to the time-ordered data to remove a low level oil/ f noise. We infer and correct 
for a small (^1%) transmission imbalance between the two sky inputs to each differential radiometer, and 
we subtract a small sidclobc correction from the 23 GHz (K band) map prior to further analysis. No other 
systematic error corrections are applied to the data. Calibration and baseline artifacts, including the re- 
sponse to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical 
uncertainties in the characterization of the beam response. Both arc accounted for in the covariance matrix 
of the window function and are propagated to uncertainties in the final power spectrum. We characterize 
the combined upper limits to residual systematic uncertainties through the pixel-pixel covariance 

This paper is organized as follows. In §1.1 we define the terms and notation used throughout the paper. 
In §2 we discuss the iterative algorithm for making maps from time-ordered data, then generalize to the case 
of simultaneous calibration and sky map estimation. (Appendix C further gcneralizc;s to map-making with 
polarization data.) We also discuss the noise properties of the time-ordered data. In §3 we discuss combined 
systematic error limits from calibration and map-making. We also present the noise properties of the sky 
maps in terms of their pixel-pixel covariance. Finally, we derive upper limits to environmental perturbations 
and summarize the combined systematic error budget. In §4 we present our conclusions. 

1.1. Notation and Overview 

Throughout this paper, we denote vectors and scalars with bold and plain lowercase letters, respectively. 
Matrices and operators are denoted with uppercase bold letters. Following Stompor et al. (2002) we denote 
vector and matrix component indices in parentheses, saving subscripts and superscripts to further identify 
quantities. A summary of the most frequently used symbols is given in Table 1. Unless otherwise stated, all 
temperatures are specified in thermodynamic units. 

WMAP measures the brightness temperature of the sky as a function of position, t{9, 0) t(p), where 
p denotes the sky map pixel number, indexed from 0, in HEALPix nested format (Gorski et al. 1998). To 
make this measurement, WMAP scans the sky and measiircs the temperatiirc difference between two points 
at time t. The resulting time-ordered differential data (TOD) is denoted d. The main goal of the map- 
making is to obtain the minimum variance estimate of the sky map, t, by inverting the raw differential data. 
Note that t is the true sky temperature convolved with the nominal instrument beam, plus instrument noise. 
In the process of solving for the map, we calibrate the data by estimating the gain and baseline from the 
fiight data itself; characterize the full instrument beam response function from observations of the planets; 
characterize the noise spectrum of the instrument, and place limits on residual systematic errors. 

In order to produce stable data with a nearly-white noise spectrum WMAP employs 20 high-electron- 
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mobility-transistor (HEMT) based differential radiometers. Each radiometer measures the brightness differ- 
ence between two inputs, one fed by an A-side beam, the other by a B-side beam approximately 141°apart. 
A detailed description of their design and fabrication may be found in Jarosik et al. (2003a); a summary of 
their in-flight performance is presented in Jarosik et al. (2003b). The 20 radiometers form 10 polarization- 
sensitive "differencing assemblies" (DA) which are designated based on their frequency or waveguide band: 
Kl, Kal, Ql, Q2, VI, V2, Wl, W2, W3, W4. The two radiometers in a DA are sensitive to orthogonal 
linear polarization modes; the radiometers are designated 1 or 2 (e.g., Kll or K12) depending on which 
polarization mode is being sensed. Each of the 20 radiometers is intrinsically a 2-channel device, with chan- 
nels designated 3 and 4 in the flight telemetry, e.g., K113 or K114. [Channels 3 and 4 were designated left 
and right, respectively, in Jarosik et al. (2003a).] There are 40 such data channels in the flight telemetry. 
As discussed below, each of the 40 channels is individually calibrated, then the 4 channels from a single 
differencing assembly are combined to form differential intensity and polarization signals as follows. Let dij 
be the calibrated differential signal from a single channel, j, of radiometer i. The differential intensity data 
(Stokes parameter /) is the average of all 4 channels 

d= ^(dl3 + di4) + ^(d23 + d24). (1) 

The differential polarization data is obtained by taking the difference between the two radiometer signals 

P = ^ (dl3 + di4) - ^ (d23 + d24) . (2) 

In Appendix C on polarization map making, we relate the differential polarization signal to the Stokes 
parameters Q and U. Kogut et al. (2003) discuss additional aspects of polarization mapping and analyze 
the flrst-year temperature-polarization correlation data based on these maps. Note that we can also form 
null channels from the data by taking channel differences, (djs — di4). As discussed in §3.2, these channel 
combinations provide valuable consistency tests for the final sky maps. 

A single channel of uncalibrated differential data may be modeled as 

c = g [M • (t + t,) + n] + b, (3) 
where each quantity is a function of time: 

c{t): raw differential data, in counts or digital units ("du"). 
g{t): instrument responsivity (here called gain), in du mK~-^. 
h{t): instrument baseline, in du. 
n{t): instrument noise, in mK. 

M • (t -h ts) = At(t): differential sky signal from all sources, convolved with the instrument beam, in 
mK. This includes fixed sources, t (e.g., CMB and Galactic emission) and moving sources, ts (e.g., 
planets). 

In practice, the differential signal is integrated for a fixed time t and sampled at discrete times ti, thus we 
may regard time series data as a vector with Nt observations. The integration time per observation is 128.0 
ms, 128.0 ms, 102.4 ms, 76.8 ms, and 51.2 ms for bands K through W, respectively. 
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The differential temperature at time t is the convolution of a time-dependent mapping function, M, 
with the sky signal at time t, t + ts{t) 

At(i) = J rfO„/ M(n(t), n') (t(n') + t,(n', t)). (4) 

Here t = tc + tg is the fixed sky signal consisting of CMB anisotropy, tc, and Galactic foreground signals, tg, 
while ts represents all time dependent sources, especially the Sun, Earth, and Moon, which are potentially 
visible in the far sidclobcs of the instrument. The operator M can be represented as an Nt x Np matrix where 
each row is the normalized, full-sky beam response in sky-fixed coordinates as given by the scan pattern. 
Several features of the mapping function that pertain to the data processing are discussed in Appendix A. 
The main beam response is mapped in flight using observations of Jupiter as a far-field point source (Page 
ct al. 2003a). An important aspect of the WMAP optical design (Page ct al. 2003b) was to limit sidelobe 
pickup to negligible levels and to have the effective beam response in the final sky maps be approximately 
circularly symmetric. We discuss each of these topics in more detail in separate papers, (Barnes et al. 2003; 
Page et al. 2003a), while this paper summarizes the main results in terms of systematic error limits in the 
sky maps. 

The instrument gain, baseline, and noise are determined from the flight data itself. This is an iterative 
process that we discuss in detail in §2.2.1. Here we provide a brief overview of our terminology in order to 
frame the following discussion of systematic errors. Let the gain and baseline measured in flight be g and 
b, respectively. The calibrated differential signal is then 

-(c — b)g gb — b 

g g g g 
With calibrated data available, the sky map is obtained by evaluating the linear equation 

t = Wd, (6) 

where W is a linear operator defined in §2.1. The properties of W are determined by the scan strategy of 
the observatory and the noise properties of the time-ordered data, d. 

The WMAP scan pattern is an integral part of the mission design (Bennett et al. 2003a). It consists of 
a compound spin and precession centered about the Sun- WMAP line, with parameters as given in Table 2. 

There are several aspects of this scan strategy that are important for high quality data: scans of a given pixel 
cross at many angles so that the effective beam response is symmetric; a given pixel is observed on many 
different time scales from minutes to months; the angular velocity of a given line of sight is nearly constant 
on the sky; the instrument observes a large fraction (>30%) of the sky each day; and the time-average of the 
differential data is approximately zero over an hourly calibration period, allowing for robust initial baseline 
estimation. 



2. THE MAP-MAKING PIPELINE 

A graphical overview of the WMAP data processing and analysis pipeline is shown in Figure 1 of Bennett 
et al. (2003b). The heart of the pipeline is a set of programs that bring science and housekeeping data from 
the Science and Mission Operations Center (SMOC) through to a set of calibrated full sky maps for each of 
the 10 WMAP differencing assemblies. Numerous additional programs are used to generate ancillary data 
products such as beam response maps, calibration files and analysis products. In this section we provide a 
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high-level overview of each program in the map-making pipeline and provide references to later sections of 

this paper, or to companion papers, for further details. 

Raw telemetry data from the satellite is transferred approximately once per day through NASA's Deep 
Space Network (DSN), to the SMOC, located at the Goddard Space Flight Center. The SMOC monitors the 
basic health and safety of the Observatory, sends all commands, and requests re-transmissions of data that 
were missing from a previous transmission. The data are then "level-0" processed into a set of time-ordered, 
daily files which contain science data, instrument housekeeping data, spacecraft data (including attitude and 
ephemeris data), and event message files. These files are then transferred to the Science Team's processing 
facility, also at the Goddard Space Flight Center. 

Every time a new full day of data arrives, a series of automated procedures perform the following tasks: 

1) Generate a standard set of daily plots that are archived and visually inspected; 2) Generate a reduced 
"trending archive" . which consists of subsets of the data sampled once every 10 minutes. In the case of the 
science data, we record the mean and rms of each channel in a 10 minute interval, whereas the housekeeping 
data are sub-sampled, once per 10 minutes; 3) Perform a series of data quality checks that search the data 
for violation of prc-sct range limits or excessive time- gradients in the telemetry signals. Limit violations are 
logged and notification is sent to a member of the science team via e-mail. (Initial limit tests are performed 
at the SMOC as well.) 

At selected time intervals, the level-0 telemetry files are collated by a pre-preprocessor into a master 
archive of raw (uncalibrated) data. The major functions of the pre-processor are to: 1) collate the science 
and housekeeping data into single daily files; 2) flag data that is Huspcx;tcxl or known to be unusable; 3) 
interpolate the attitude and ephemeris data to times that are commensurate with the science data time 
stamps; and 4) apply a coarse flag to data that is within 7° of one of the outer planets (Mars through 
Neptune) so that it may be rejected from the initial sky maps, but identified for beam mapping (Page et al. 
2003a). (Only Jupiter data is used for making the final beam maps.) 

Initial sky maps and calibration data arc generated from the raw archive using the iterative map-making 
algorithm first described in Wright et al. (1996), and further described in §2.1. As discussed in §2.2, the initial 
calibration is determined by fitting the raw time-ordered data to the known signal produced by the CMB 
dipole. Because the sky signal contains significant higher-order power (Z > 1), the calibration solution must 
be iteratively improved in concert with the initial sky map iterations. The convergence of this simultaneous 
fit has been demonstrated with end-to-end simulations, which are also described in §2.2. 

The calibration solution converges more rapidly than the sky map does, so we freeze the initial calibration 
solution after ^10 iterations before proceeding to convergence with the sky maps. With the initial dipole- 
based calibration data in place, a re-processor generates a refined gain and baseline solution and applies this 

to the data. The program also updates the data quality flags, as necessary, then writes the calibrated data to 
a new final time-ordered archive. The refinements to the dipole-based gain solution are discussed in Jarosik 
et al. (2003b) and in §2.2.1. The initial baseline solution is refined with a pre-whitening filter (Wright 1996) 
which is presented in detail in §2.3.2. 

The final sky maps are computed using the final calibrated time-ordered data as input. The first-year 

sky maps required 20 post-calibration iterations to be sufficiently converged. The map-making algorithm is 
fundamentally the same as is used in the initial estimates, but we add some refinements for this final stage 
of processing. These include: 1) correcting for a small (Sl%) loss imbalance between the A and B-side 
sky beams. Jarosik et al. (2003b) demonstrate that this effect is nearly orthogonal to the gain solution, so 
its inclusion after the calibration processing does not invalidate the gain solution. 2) We weight individual 
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time-ordered observations by their proper statistical weight to account for the small change in instrument 
noise (<1%) over the course of a year due to the 0.9 K seasonal temperature variation experienced by 
the instrument cold stage. 3) We compute the planet-boresight angle for each observation to minimize 
unnecessarily conservative data loss. The criterion used for the first-year maps is a cut of radius 1?5. 

We also generate polarization maps using a generalization of the temperature algorithm, the details 
of which are presented in Appendix C. For the first-year data processing, we have generated maps of the 
Stokes parameters Q and U , but we have not yet fully characterized all of the potential systematic errors in 
these maps. However, the temperature-polarization correlation data are much less prone to systematic errors 
than the polarization auto-correlation data. Kogut et al. (2003) have analyzed the temperature-polarization 
data and, supported by systematic error limits from Barnes et al. (2003), they find a significant correlation, 
including the signature of a relatively early epoch of cosmic reionization. 

In the remainder of this paper, we present a more detailed description of each of the map-making and 
calibration procedures, including an assessment of their performance with the first-year WMAP data. We 
then derive detailed systematic error limits applicable to the time-ordered data, to the sky maps, and to the 
angular power spectrum. 



2.1. Map-Making with Pre-Calibrated Data 

While the process of generating the final sky maps from calibrated data comes last in the map-making 
pipeline, we discuss the algorithm first because the algebra of map-making is central to the entire data 
processing scheme, and it helps to guide the systematic error analysis. 

We consider the problem of estimating a sky map, t, from calibrated, differential time-ordered data, d, 
which is a linear function of the sky map 

d = Mt + n, (7) 

where M is the mapping function of the experiment, which has Np columns and Nt rows, Np is the number 
of sky map pixels and Nt is the total number of time-ordered observations. In its simplest form, each row 
(observation) of the scan matrix contains a in the column (pixel) seen by the A-sidc beam, and a —1 in 
the column (pixel) seen by the B-side beam. This matrix can be generalized to include the effects of beam 
convolution, but for WMAP these refinements are small and are being deferred to future processing. The 
effects of a differential loss imbalance between the A and B-side beams is readily accounted for by using 
values different from ±1 in M. An analysis of this effect in the WMAP radiometers is presented in Jarosik 
et al. (2003b). The details of how we account for it in the pipeline are given in Appendices A and C. 

The noise n is assumed to have zero mean and covariance N, 

(n) = (8) 
(nn"^) = N. (9) 

We defer a detailed discussion of the WMAP noise properties to §2.3.2, but for most radiometers it is rea- 
sonable to approximate the noise covariance as diagonal, N ~ (Tq I, where ctq is the rms noise per observation 
and I is the identity matrix, though this assumption is not required for the algorithm described below to 
converge. 

The least-squares, or maximum-likelihood estimate of the sky map, t, results from solving the normal 
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equations 

t = (M^N^iM)"^ • (M^N^M). (10) 

More generally, we obtain an unbiased estimate of the sky map by choosing any symmetric matrix S in place 
of (Tegmark 1997). To see this, substitute Mt + n in place of d in the above equation to get 

t = (M^SM)-i • (M^S [Mt + n]) = t + (M^SM)-i • (M^Sn). (11) 

Thus t reduces to t plus a noise term that is independent of t and has zero mean over an ensemble average. 
For the first-year processing we take S = I. To simplify notation, we define a matrix W = (M-^M)"^ • 
in which case the map solution is t = Wd. 

The pixel-pixel noise covariance in the sky map solution is 

S = ((t-t)(t-t)'^) = ((Wn)(Wn)'^) =WNW'^ 

= (M'^M)-^ • (M'^NM) • (M'^M)-^ (12) 

In the limit that N is diagonal and the rms noise per observation, uo, is constant, S reduces to (M^M)~^. 

For the WMAP scan pattern, the matrix M-^M is diagonally dominant with diagonal elements Uohsip), the 
number of observations of pixel p by either the A or B-side beam. Thus, to a very good approximation, the 
pixel-pixel covariance matrix is diagonal 

2 

nobs(Pi) 

Values for ctq are given by Bennett et al. (2003b). 

The leading order off-diagonal terms occur at the beam separation angle (^beam ^ 141°), and are of 
order 0.3% of the diagonal elements. If the time-ordered noise N is not diagonal, then maps produced with 
the above algorithm will have correlated noise (stripes) along the scan paths defined by M. This is a small, 
but not negligible, effect for some of the WMAP radiometers, and is entirely negligible for others. The noise 
properties of the time-ordered data and sky maps are further discussed in §2.3.2 and §3.2, respectively. 



2.1.1. Iterative Map Making 

The evaluation of the sky map solution Wd requires the inversion of the Np x Np matrix D = M^M. 
We use the iterative approach introduced by Wright et al. (1996) to evaluate this expression. Briefly, suppose 
we have an initial guess for the sky map, to, which differs from the true sky map, t, by Sto = to — t. Then 
D to = D (t -h (5to) can be recast as 

D Jto = D to - M'^d, (14) 

where we have used the fact that Dt = M'^d. As noted above, D is diagonally dominant, so a good 
approximate inverse for D is 

Ilobs \Pi ) 

This leads to the approximate solution for the residual, dt^ ~ D^^[Dto — M^d], and suggests the following 
iterative solution 

tn-l-l = t„ - 6tn (16) 

= t„-D-i[Dt„-M^d] (17) 
= (D-iM^)d+(I-D-iD)t„. (18) 
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The interpretation of equation (18) is that for each pixel the new sky map temperature is the average of all 
differential observations of that pixel (accounting for the sign of the observing beam) corrected by an estimate 
of the signal in the paired beam, based on the previous sky map iteration. The expression in equation (18) 
can be efficiently evaluated because the sums can be accumulated by reading through the time-ordered data 
from disk, each iteration, and accumulating data into arrays of length Np 



where = 1 in the initial sky map processing and is proportional to a noise weight (equation C31) in the 
final sky map processing. Note that it is never necessary to store or invert an Np x Np matrix. 

We have tested this algorithm extensively with fiight-like simulations. In this section we present results 
for an "ideal" noiseless instrument with circular beams and perfect calibration to isolate the performance of 
the map-making algorithm from other effects. More realistic data models are introduced to the simulation in 
subsequent sections. Figure 1 shows a sample residual map, tout — Un, generated from a one- year simulation 
of Q2 data. The input sky map included realistic CMB signal with a peak-to-peak amplitude of ^ ±420 
/xK, and a Galactic signal with a peak brightness of ^50 mK. The output sky map is recovered with an 
rms error of < 0.2 /xK, after 50 iterations. The dominant structure in the residual map is a mode aligned 
with the ecliptic plane. The power in this mode is concentrated in spherical harmonic mode I = 4, due 
to a combination of the WMAP scan strategy and the beam separation angle. This is the mode on the 
sky that is least well measured by WMAP (except for the monopole!) and is thus the slowest to converge, 
though additional iterations would reduce its amplitude even further. The final first-year fiight sky maps 
were effectively run for 80 iterations. Since the rms error associated with this term is very small, and since we 
build up a more realistic data model in subsequent simulations, we do not further quantify this contribution 
to the final systematic error budget. Rather, we subsume it into an overall map-making and calibration error 
budget that includes this and other effects together. 

This iterative approach to map-making is readily generalized to polarization maps as well - the formalism 
is presented in Appendix C. We have tested that algorithm with the same simulations used above to test 
the temperature algorithm and find that the polarization maps converge even faster than the temperature 
maps. After 10 iterations, the map-making artifacts in a residual polarization map are <0.05 fxK peak-peak. 



The processing algorithm described above assumes that the data have already been calibrated. In 
practice, wc; use the above; algorithm in the second stage of map-making, after an initial stage in which we 
simultaneously solve for the radiometer calibration and the sky map. In the initial stage of map-making, 
we employ the same iterative algorithm to solve for the map, but rather than processing straight through 
the time-ordered data on each iteration, we process the data one hour at a time, pausing to solve for the 
calibration in each radiometer channel, before accumulating the calibrated data. The calibration solution 
then iteratively improves as the sky model improves. The following sub-sections lay out the procedure in 
detail, and present results for the flight data with an assessment of its precision and accuracy based on 
flight-like simulations. 



"obs • tn+l{PA) Hobs • t„+i(pA) + Wi [d{ti) + t„(ps)] 
Hobs • tn+liPs) Hobs • tn+liPs) - Wi [d(ij) - t„(pA)] 



nobs(PA) nobs(PA) + Wi 

"obsCPs) nobs(ps) + Wi 



(19) 



2.2. Combined Calibration and Map Making 
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2.2.1. Instrument Calibration from the Dipole Modulation 

For a sufficiently short period of time the instrument gain and basehne can be approximated as constant, 
c ^ .9fe(At + n) + bk, where gk and bk arc the gain and basehne during the /c**^ 1-hr cahbration period. Since 
the sky signal At is dominated by the CMB dipole measured by COBE, Atd, a single channel of raw data 
can be modeled as 

Cm(5fe,&fe) =5'fc(Atd + Atv) + 6fe, (20) 

where Atv is the additional dipole moment induced by the motion of WMAP relative to the solar system 
barycenter (the rest frame of the COBE dipole). 

We fit for the gain and baseline in each calibration period, k, by minimizing 



^2 ^ ^ Htj) - C^{t^\gk,bk)] ^^^^ 

where i is a time-ordered datum index. We omit data that are flagged as unusable, and data when either 
the A or B-side beam points within a Galactic pixel mask. The mask used for this latter application is the 
Kp8 mask defined in Bennett ct al. (2003c), without edge smoothing. This mask is used throughout the 
map-making pipeline. The fit is performed for each of the 40 WMAP channels independently. To minimize 
the covariance between the recovered gain and the baseline, it is necessary to have a scan strategy such 
that the time average of the sky signal, At, is approximately zero in one calibration period. The combined 
spin and precession of WMAP is designed to produce time-ordered data that satisfies this requirement. For 
example, in K band, which has the largest sky signal, a 1-hr running mean of the differential sky signal has 
an rms fluctuation of 14 ^K, compared to a dipole signal of greater than 3 mK. After each hour of data is 
processed for the calibration solution, the data are accumulated as per equation (19) to develop the sky map 
solution. 

The largest source of error in the calibration fit is due to un- modeled sky signal from the CMB anisotropy 
and Galactic foreground emission, Atg = At — Atj. This projects onto the dipole signal and, as shown 
below, causes errors in the gain solution as large as 5-10% in K band, where the Galactic signal is strongest. 
The calibration fit may be iteratively improved by subtracting an estimate of the anisotropy from the raw 
data prior to fitting. In particular, let g[, be the gain inferred for calibration period k from the previous 
iteration of the calibration fit. Then minimize 

^2 ^ {c'jti) - Clnit^\gk,bk)f ^22) 



2 



where 



ak^t'. (23) 



and Atg is the differential sky signal (less the dipole component) computed from the previous sky map 
iteration. This process is repeated until the calibration solution is sufficiently converged. 

Note that the absolute calibration is tied to the time-dependent portion of the dipole signal, Atv; we 
use the fixed dipole as a short-term transfer standard only. In particular, when we update the sky model 
and apply the anisotropy correction in equation (23), any error in the fixed dipole moment, Atd, is assumed 
to be anisotropy, and is applied as a correction in the same way. For a data set of at least one year in length 
(one full cycle of Atv), the error in the absolute calibration will be essentially orthogonal to any error in 
the fixed dipole Ata- 
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2.2.2. Performance of the Dipole-Based Gain Solution 

As an illustration of the systematic gain error induced by higher-order (l > 1) anisotropy, Figure 2 
shows an example of the gain solution convergence from a one-year low-noise simulation. This simulation 
implements the simultaneous calibration and sky map estimation discussed above and was run for 30 iter- 
ations. The example shown is for one channel of K band data (the worst case) which exhibits ~7% errors 
after one iteration, corresponding to a sky model that has only a dipole component. After 30 iterations, 
the residual errors are <0.1% over the entire year. Similar, or better, performance is achieved for all other 
WMAP channels. 

In processing the final first-year maps, the combined calibration and map-making code was run for 10 
iterations. However the initial sky model was based on an earlier "pathfinder" run of the pipeline that 
ran for a total of 30 iterations of combined calibration and map-making plus an additional 20 iterations of 
sky map convergence. Thus we conservatively estimate that the combined absolute and relative calibration 
errors due to incomplete calibration convergence to be <0.1%. We defer a discussion of the final calibration 
uncertainty to §3.1. 

Figure 3 shows a sample of the converged gain solution from equation (22) for two WMAP channels, 
K113 and V113. Note that the radiometer gains are typically drifting by a few percent over the course 
of the first year. As wc show below, the dipolc-bascd fits easily track drifts at this level. The noise in 
the gain solution is typically a few percent per hourly calibration period, though, as is readily seen in the 
figure, the noise level changes with time of year as the scan pattern sweeps around the CMB dipole. The 
VI 13 gain exhibits an additional modulation that is clearly correlated with the physical temperature of the 
instrument. However, the time scale of the temperature change is slow enough that the corresponding gain 
changes are well tracked by the dipole fits. Quantitative limits on thermally induced gain and baseline errors 
are discussed in Jarosik et al. (2003b) and in §3.4.1. A summary of the gain statistics from the flight data 
is given in Table 3. 

2.2.3. The Initial Baseline Solution 

The bottom two panels of Figure 3 show the converged baseline solution resulting from the fit in 
equation (22) for one year of K113 and VI 13 flight data. The fits have had a mean subtracted, and have 
been divided by the gain to convert to temperature. These plots, which are representative of all 40 channels, 
show that the offsets of the radiometers arc typically stable to ±5 mK over the course of the first year. 
Simulations demonstrate that the hourly baseline solution is unbiased. However, it is also clearly noisier 
than optimal, consistent with the flight measurements of the noise power spectral density (Jarosik et al. 
2003b). In §2.3.2 we describe an improved baseline model that is based on the application of a pre- whitening 
filter tailored to the measured noise spectrum of each channel. 

2.3. Improving the Calibration Model 

The sky maps obtained with the hourly calibration are reasonable; however the noise in the calibration 
solution, particularly in the baseline, is significantly higher than optimal, and the use a piecewise continuous 

calibration in the final maps would introduce striping in its own right. In the following subsections, we 
present the steps undertaken to filter the gain and baseline solutions that enter into the final sky maps. 



Prior to generating the final maps, this refined caUbration is appUed to the data and written to disk as the 
final first year calibrated time-ordered archive. 



2.3.1. The Gain Model 

Jarosik et al. (2003b) present a physical model for the gain that is based on the RF bias, or "total 
power" measured in each channel, and on the physical temperature of the instrument cold stage, which is 
monitored with high resolution platinum resistance thermistors (PRTs) . Each of these quantities is recorded 
once every 23 s in the engineering telemetry with a relative noise that is substantially lower than the noise in 
the dipole-based gain solution. Thus, if the model fits the dipole-based data satisfactorily, it offers a means 
for measuring the gain with more precision, and on time scales shorter than the spin period. The model for 
the gain, g(t), has the form 

where V is the measured RF detector bias, TppA is the measured temperature of the FPA, and go, Vq, and 
To are fit constants. See Jarosik et al. (2003b) for more detail. 

Figure 8 in Jarosik ct al. (2003b) shows the performance of the gain model when fit to the dipole-based 
gain solution (see also Figure 7 in this paper). In §3.1 we evaluate the overall performance of the gain model 
and the hourly gain solution in the context of an end-end simulation designed to place limits on combined 
calibration and map-making errors. We will conclude that the model provides an excellent description of 
the radiometer gain, and here we adopt it as the final gain solution for further processing. The gain model 
fits into the data processing sequence as follows. After we iterate the simultaneous calibration and sky map 
solution long enough for the calibration to converge (10 iterations when starting with a good sky model), we 
freeze the dipole-based calibration and fit the gain model parameters in equation (24). All subsequent data 
products are produced with data calibrated using this gain solution, including the time-ordered archive, the 
final sky maps, and the Jupiter beam maps. 



2.3.2. Baseline Filtering 

The baseline that results from the initial calibration solution is not optimal. This is due to the fact 
that the initial baseline is sampled once per hour (0.28 mHz), while Jarosik et al. (2003b) show that the 
power spectral density of the noise has a 1// knee frequency of a few mHz, typically. If the initial baseline 
estimate were used in the final sky maps, it would generate weak stripes of correlated noise along the scan 
paths, as per equation (12). Even so, it is important to note that 1// effects are small relative to the white 
noise. In the worst WMAP radiometer, W41, the amplitude of the noise covariance N at small lag is ^2% 
of the white noise variance. Thus we treat 1// noise iteratively in the data processing by first ignoring it to 
obtain an estimate of the gain, baseline, and sky solution. Then we subtract the estimated sky signal from 
the time-ordered data, apply a pre-whitening baseline filter to the residual noise, add the sky signal back 
in, and write the data to a final calibrated, time-ordered data archive. The approach of first subtracting an 
estimated sky signal is designed to avoid biasing the gain solution and/or removing low-order power from 
the sky maps. The noise properties of maps constructed in this way must account for the filtering process. 
We discuss the map-making algebra appropriate to our filter implementation in Appendix B. 

The steps we follow to define and apply the filter are as follows. We remove an estimate of the sky 
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signal, in du, from the raw differential data using 



c'{t) = c{t) - g(t) At'(i) - bk, 



(25) 



where bk is the hourly baseline point appropriate to the current time, g is the final gain solution from 
equation (24) and At' is the differential sky signal computed from the initial sky map. We then evaluate 
the auto-correlation function of c'{t) to a lag of 10^ sec. Results for representative radiometers are shown 
in Figure 4. This range of lags is sufficient to account for both the long-range correlations due to 1// noise 
and the correlation at a lag of 1 observation due to the low-pass post-demodulation filter in the Analog 
Electronics Unit (Bennett et al. 2003a). The baseline filters are then defined as follows (Wright 1996) 

1. Fit the auto-corrclation function, C(At), to the model defined below. 

2. Fourier transform the model correlation function to generate the power spectral density P{f). 

3. Compute the Fourier space filter w{f) = l/^/P{f) and set w{0) = to produce a zero mean output 



4. Fourier transform the filter w{f) to generate the time domain filter w{t), normalized to 1 at t = 0. 
By inspection, the auto-correlation functions are well modeled by a log-linear function 



where r is the integration time for a single observation, Ci is the correlation at lag r, measured from the 
data, and A and B are fit parameters. Note that A gives the typical fractional covariance at small lag, while 
the suppression of correlations at large lag (~2000 s) is dictated by the subtraction of the hourly baseline 
as a pre-filter. The best-fit parameters are given in Table 4, and fits for selected radiometers are shown in 
Figure 4. 

The derived pre-whitening filters, w{f), are plotted as a function of frequency in Figure 4. One point 
of particular interest is the filter response at the spin frequency, 7.7 mHz. As shown in Table 4, the best 
channels have a transmission of > 95%, while the worst case, W41, is just above 35%. These values indicate 
the amount by which the dipole (calibration) signal would be suppressed if the filter were applied prior to 
calibration and sky signal subtraction. The convolution of c'{t) with 'w{t) is performed in Fourier space 
using the Numerical Recipes routine convlv (Press et al. 1992). The number of data points convolved at 
any one time is chosen to be the smallest power of 2 such that the data span a full day with sufficient padding 
beyond the day to guarantee that wrap-around effects are negligible. This is 2^° for K-Q bands, and 2^^ for 
V,W bands, which gives a minimum of 2.9 hours of padding on each end of a day. Sample auto-correlation 
functions obtained from the filtered data are shown in Figure 4. The filtering is clearly effective at removing 
low frequency noise in the time-ordered data. Another example of filtered data is seen in Figure 5, which 
shows 1 day of W42 data, one of the worst radiometers for 1// noise, before and after filtering. These data 
are smoothed with a 46 s window to show structure in the unfiltered data since plots of unsmoothed data 
before and after filtering arc virtually indistinguishable. 

The above results are encouraging but not definitive, because the process of sky signal subtraction and 
re-addition could introduce correlated artifacts that these tests would miss. The ultimate test of a filter is 
its ability to "clean" the pixel-to-pixel covariance matrix of the final sky maps and the noise covariance of 



signal. 




At = T 

T<At< 10^/-^ s 
At > 10^/-^ s, 



(26) 
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the angular power spectrum, without altering the underlying sky signal. The sky map noise properties are 
discussed in §3.2, while the noise properties of the power spectrum are quantified in Hinshaw et al. (2003b). 

2.3.3. Baseline Jumps 

Limon et al. (2003) identify 21 instances of sudden baseline jumps, or "glitches" , during WMAP^s first 
year of operation. These events have been identified as small shifts in the properties of several microwave 
components resulting from sudden releases of internal mechanical stress, presumably from thermal changes. 
These events last for less than 1 s, and cause no discernible change in the radiometer gain or noise properties. 

Care must be taken in the application of the baseline filter in the vicinity of these steps to avoid ringing 
in the filtered data. Each event is initially flagged by the pre-processor for at least 1.2 hours on either side of 
the event. Since the initial hourly baseline is derived entirely in the time domain, the ±1.2 hr flagged interval 
ensures that this baseline estimator only "knows" about data on one side of the jump or the other. Prior to 
convolving the raw data with the pre-whitening filter, we subtract the initial hourly baseline from the data 
as a pre-conditioner. Thus all data that is input to the convolution routine has approximately zero mean. 
On output, the re-processor expands the flagged interval by 0.5 hr on either side of the event to ensure that 
no edge effects propagate into the usable data. In the first year of operation, a total of 0.13% of the data 
was lost to these steps. See Table 2 in Bennett et al. (2003b). 

The threshold amplitude for jump detection by visual inspection is ~0.05 du, which corresponds to a 
jump of ~150 jjK in the calibrated output of radiometer W12, the worst off'ender. To assess the effect of 
undetected baseline steps in the data, we have generated a test data set in which we take 24 hours of flight 
W12 data and insert a step of 0.05 du in each channel. We then run the data through the pre-whitening 
filter to see the magnitude of the transient response. The result is a transient baseline error with a peak 
magnitude of 80 /xK, which lasts for less than one 2-minute spin period. The total time the baseline error 
exceeds 10 /iK is 22 minutes, or approximately 11 spin periods. We pessimistically assume that there could 
be as many as 40 steps at or just below the threshold of detection, and that half of these are in W12. If 
we assume these occur at random times, and note that WMAP observes ^30% of the sky in any given 
hour, then any given sky pixel is likely to "see" approximately llx20x0.3Ri66 data points with baseline 
errors greater than 10 /iK. Since the sign of a given step is random, and since Wll data is combined with 
W12 in the sky maps, we estimate the residual systematic error in a given pixel of the Wl map is less than 
10 ^K /2/\/20 ~ 1 /uK. We emphasize that no jumps have been observed in any other W band radiometer, 
thus DA-DA consistency is an excellent test of whether any statistic is sensitive to baseline errors of this 
nature. We have found no evidence that the Wl map is "out of family" (§3.2). 

2.4. Final Sky Map Processing 

Once the calibrated, time-ordered archive has been written, final sky map processing commences based 
on the algorithm presented in §2.1. At this stage in the sky map processing, we add a few features to the 
algorithm that, for simplicity, are not present in the combined calibration and map-making code. These 
include: 1) Weighting each datum by a true weight, l/u^ based on an estimate derived from the physical 
temperature of the instrument cold stage. This introduces ~1% variations in the data weights over the year, 
since the instrument noise is a weak function of temperature, and the temperature varies by ~1% over the 
year. 2) Accounting for loss imbalance, as discussed in Appendix A. In effect, we model each differential 
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observation as At = {l + Xirn)t{pA) — (^ — Xim)t{pB), where Xim is the smah loss imbalance parameter given 
by Jarosik et al. (2003b). 3) Computing the planet avoidance flag at run time to reduce the amount of data 
lost. In the final sky maps, a total of 0.11% of the data was lost to planet avoidance. See Table 2 in Bennett 
et al. (2003b). 

The final stage of sky map processing, based on the filtered data, was run for a total of 20 iterations. 
Convergence was determined by measuring the rms difi^erence between pairs of iterations for a given differ- 
encing assembly. For example, the difference between the 10th and 20th iteration of the W2 sky map is 0.08 
nK rms. We estimate that artifacts due to lack of solution convergence are <0.1 iiK rms with all of the 
power being in the low multipoles, I < 10. We present a final combined estimate of sky map artifacts due to 
calibration and map-making errors in §3.1. This estimate includes the convergence limits given above. 

To assess the effect of the improved calibration model on the final sky maps, we form differences between 
the final post-filtered maps and the last iteration of the intial, prc-filtcrcd maps. The results for DAs W3 
and W4 are shown in Figure 6. The top panels in this Figure show the difference maps from a one-year 
simulation (§3.1) that included a realistic radiometer noise and gain model. The bottom two panels show 
differences from the flight data. Because these maps are largely based on the same data, most of the white 
noise drops out of these differences. The remaining "blobs" of white noise result from the change in the 
planet cut and can be ignored. The striking feature is the striping present in the W4 difference, but virtually 
absent in the W3 difference. As we show in Appendix B, this is the structure that has been removed from 
the data by the pre-whitening filter, an interpretation that is substantiated by the analysis of the simulation. 
The fact that the W3 difference is very small is an indication that the level of striping in the unfiltered W3 
data was very small to start with. We estimate the level of residual striping in the final maps in §3.2. Images 
of the final maps at each frequency are presented by Bennett et al. (2003b). 

3. SYSTEMATIC ERROR ANALYSIS 

As discussed in the introduction, systematic errors may be classified by the nature of their source. In 
this section we place limits on the level of systematic errors in the final sky maps, using that classification 
to guide the analysis. In §3.1 we place limits on combined calibration and map-making artifacts, based 
largely on a detailed simulation of the first year of WMAP operation. In §3.2 we present null tests based on 
difference maps formed from a variety of data combinations, each of which should yield no sky signal. We 
use these maps to measure or place limits on correlated pixel noise (striping) in the final first-year maps. 
In §3.3 we discuss systematic errors relating to pointing and beam mapping errors. We conclude by placing 
stringent limits on residual errors due to environmental (thermal and electrical) and other miscellaneous 
sources. 



3.1. Calibration and Map-Making Errors 

To assess the combined errors from calibration and map-making artifacts, we have generated a high 
fidelity simulation that includes all of the effects we believe are important for calibration and map-making. 
In particular this simulation includes: 1) A sky model that closely mimics the statistical properties of the 
observed sky; 2) A realistic noise model for every channel, including 1 // noise [see Jarosik et al. (2003b) for 
a tabulation of 1// knee frequencies]; 3) A model for the thermal drift of the gain, baseline and offset of 
each radiometer, based on measured susceptibility coefficients, and driven by the actual temperature profile 
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measured in flight. This simulation generates the sky signal using a circular beam approximation. The 
effects of elliptical beams are treated in a separate, noiseless simulation in §3.3.4. We write simulated science 
data to files that mimic the raw telemetry, then process the data using the same pipeline as was used to 
process the flight data. 

The top panel of Figure 7 shows the converged gain solution from the simulation for channel V113; the 
bottom panel shows the corresponding result from the first year of flight data. In both panels, the "noisy" 
black traces arc the hourly gain data, binned in 24-hr samples to reduce the noise, and the green traces are 
the best-fit gain model (§2.3.1). For the simulation, the input gain used to generate the data is shown in 
grey. The absolute gain is recovered in the simulation to better than 0.1% for all 40 channels. 

The dipole signal seen by an observer moving with speed v relative to the rest frame of the CMB is 
To v/c, where Tq is the absolute temperature of the CMB, and c is the speed of light. Thus, additional sources 
of error that could affect the absolute calibration of the WMAP data include errors in the determination of 
WMAP's velocity with respect to the solar system barycenter (the point of reference for the COBE dipole) 
and errors in the absolute temperature of the CMB. The velocity of WMAP is routinely measured with 
respect to geocentric incrtial coordinates (GCI) with an accuracy of <1 cm s~^. The velocity of the Earth is 
determined from the JPL ephemeris with similar accuracy. The combined uncertainty from velocity errors is 
0.1 nK. The uncertainty in the absolute temperature of the CMB is 0.1% (Mather et al. 1999). Combining 
these uncertainties with the results of the simulation, we conservatively estimate an absolute calibration 
error of 0.5% for the first-year WMAP data. 

Errors in relative calibration can produce structure in the sky maps, beyond an overall normalization 
factor. The largest relative discrepancy between the dipole gain solution and the gain model in the fiight 
data is ^0.4% in K band, and ^0.2% in the other bands. Similar deviations are seen in the simulation, thus 
we use the simulation as our primary tool for placing systematic error limits due to relative calibration and 
map-making errors. We have generated residual maps from the simulated data by subtracting the known 
input sky signal from the maps produced by the pipeline. These residual maps exhibit no visible structure 
aside from the pixel noise. In order to assess the errors due to map-making artifacts, we compute the angular 
power spectrum, C;, of the residual maps and search for features in the spectra beyond a simple flat, white 
noise spectrum. The residual spectra for all 10 differencing assemblies are shown in Figure 8 and summarized 
in Table 5. In general, the spectra are consistent with white noise over a wide range of multipole moments, 
but clearly show the most variation at low I. Because of this, we specifically highlight these modes in Table 5, 
where we give C2, (Cj)3_io) and (C()ii_ioo fo^' ^ach of the DA's. For combined systematic error Hmits due 
to calibration and map-making, we assign twice the excess variance in each I range relative to the white 
noise plateau, cr*^* = 2 |(C';)i3a„d — (^i) 700-1000 1 • These values are also quoted in Table 5. For comparison, 
the average power in the CMB in each band is C2 ^ 130 /uK^, {Ci)^_-^q ~ 150 /uK^, and (Ci)ii_ioo ~ ^ /uK^. 

Because the simulation includes realistic models oil/ f noise and long-term thermal effects, these limits 

also implicitly limit artifacts at low I due to these effects. As we demonstrate in subsequent section, we 
feel this simulation captures all of the important radiometric characteristics of the instrument. Potential 
artifacts due to optical effects, especially pickup through the far sidelobes, are treated in Barnes et al. (2003), 
and are summarized in §3.3.5. 
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3.2. Difference Maps and Noise Properties 

Difference maps are combinations of the data that, ideally, should contain no sky signal. They provide 
insight to potential systematic errors and can be used to characterize the noise properties of the sky maps. 
The first set of difference maps we generate are between DA pairs with the same frequency and beam response, 
namely i(Ql-Q2), i(Vl-V2), and i(W12-W34), where W12 = i(Wl+W2), and W34 = i(W3+W4). 
Images of these difference maps are shown in Figure 9, along with low resolution versions of the sum (signal) 
maps to give a sense of the signal strength in each map. Aside from the pattern of the noise, which follows 
the sky coverage [see Figure 3 of Bennett et al. (2003b)] , the only visible structure in these difference maps 
is in the Galactic plane, especially in V band. This is understood to be a result of a small difference in the 
effective center frequency of the VI and V2 differencing assemblies (Jarosik et al. 2003a). In particular, the 
VI map has an effective frequency approximately 1 GHz lower than V2. Since the spectrum of the Galaxy at 
V band follows Ia(j^) ~ (Bennett et al. 2003c), we expect the Galactic signal to be ~3% brighter in VI 
than V2, which is consistent with the residual signal seen in the difference map. [A complete tabulation of 
effective center frequencies, radiometer by radiometer, is given by Jarosik et al. (2003a) for diffuse sources, 
and by Page et al. (2003a) for point sources.] Note that because the data are calibrated using the CMB 
dipole, there should be no residual CMB signal in such a difference map. A more sensitive comparison of the 
single DA maps is afforded by comparing their angular power spectra. In that case, it is easier to compare 
ac;ross frequencies because differences in beam response are readily accounted for by deconvolution. See 
Hinshaw et al. (2003a) for such a comparison. 

We generate three additional sets of difference maps using different combinations of the 4 channels 
within a DA. Specifically, we form the differences in the time-ordered data then generate maps as follows 

-(di3 + di4) - -(d23 + Ci24) Ai2, 
-(di3 - di4) + -(d23 - d24) A34, 

^(di3 -di4) - ^(d23 -d24) ^ A1234, (27) 

where indicates the map-making process. The A12 maps are based on the polarization data, but pro- 
cessed as temperature maps, i.e., without attempting to demodulate the polarization signal. Since the two 
radiometers within a DA have completely independent detection chains, and since the polarization signal 
is weak to begin with (and is further suppressed by the lack of demodulation) the noise properties of the 
A12 maps should be virtually identical to the nominal signal maps. The A34 and A1234 maps are based 
on channel differences, (d^a — d^), and since the two channels within a radiometer have partially correlated 
noise, the noise properties of these latter maps will be different than the maps based on (di3 -|-dj4). However, 
these maps do provide a check on the channel calibration, common-mode thermal effects and other potential 
artifacts. 

For each difference map the two-point correlation function and the angular power spectrum are calcu- 
lated. The results are shown in Figures 10 and 11 and summarized in Table 6. Figure 10 shows the two-point 
function computed from the A12 maps for Q2, V2, and W2. The most apparent feature in each of these 
functions is the slight bump at the beam separation angle of ^beam ~ 141°, as expected (§2.1); the first data 
column of Table 6 gives C(0boam)/C'(O), for each DA. Note that, with the exception of K band, the ratio 
is typically 0.3%. The larger K band values arise because Galactic leakage in these difference maps is most 
severe in this band. This is also the source of the weak large-scale feature in the V2 two-point function in 
Figure 10. While this residual signal is small compared to the temperature signal, it is a systematic error 
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that must be accounted for in the analysis of polarization data (Kogut et al. 2003). 

Figure 11 shows the angular power spectra of the difference maps for each of the 10 DAs, as well as for 
the final signal maps. Table 6 summarizes their statistics. Note that the white noise plateau in the 4 channel 
combinations per DA divide into two families, as noted above, due to the correlations between channels 
3 and 4. As a result, the null combinations, A34 and A1234, cannot be used to estimate the white noise 
parameter ctq for the signal maps. However, the polarization channel is seen to be in excellent agreement 
with the temperature channel in the white noise tail, thus, to the extent that real polarization signals, or 
other systematics, such as bandpass mismatch are not important, these maps should provide an excellent 
noise model for the temperature data. Table 10 summarizes the shape of the angular power spectrum at 
low I in the same way Table 5 did for the simulation. We find the spectra of these difference maps to be 
remarkably flat, with residual quadrupole moments of <4 /xK^ for all bands except K (in which the difference 
is dominated by bandpass mismatch) and a single combination of W3. This value is much smaller than the 
small quadrupole measured in our sky (Bennett et al. 2003b). The power in the other multipole ranges is 
very close to the white noise floor, as seen in the final columns of Table 6. Since the residual signals seen 
in the flight difference maps are somewhat lower than those seen in the simulation (Table 5), we adopt the 
more conservative limits from the simulation as systematic error limits for structure at low I. This allows 
for the possibility than some of the error seen in the simulation comes from, e.g., common-mode calibration 
errors that cancel in the difference maps. 

The two-point correlation function of the A 12 maps demonstrates that the angle- averaged off-diagonal 
terms of the pixel-pixel covariance matrix are less than ~0.3%. However, the maps in Figure 6 indicate 
the potential for stripes along the scan directions for which the covariance can be locally higher than the 
angle- averaged value. In order to determine the magnitude of the covariance along the scan directions we 
perform the following computation. We form W band difference maps: Wi—Wi', where i = 1-4, and Wi' is 
the average of the 3 other W band maps, e.g., Wl' = i(W2-hW3-hW4). We then form time-ordered data 
from this map using the pointing appropriate to DA Wi and compute the auto-correlation function, C{At), 
from 30 days of data. This provides a measure of the pixel-pixel covariance along a stripe. The results for 
W3 and W4 are shown in Figure 12; the top panels show the covariance, normalized to one at lag zero, 
computed from the unfiltered maps, while the bottom panels show the results for the filtered maps. In order 
to convert the time axis to angular displacement along a scan, recall that the 2?784 s~^spin rate translates 
to a 2?6 s~"'^rate for the beams in either focal plane (the second decimal place depends on position in the 
focal plane, and time in the precession cycle). The W3 result shows a hint of covariance (~0.1%) at lag 0.1 s, 
or 0?26, prior to filtering, but none after (<0.05%). Prior to filtering, the W4 result shows clear covariance 
of up to 0.5% at small lag, decaying to <0.1% at lags of ^ 10^ s, roughly one full spin. After filtering, 
the covariance is reduced by nearly a factor of two, but is still clearly detectable. This is understood to be 
residual covariance that survives the filtering process because of the fact that we subtract an estimated sky 
signal, based on the data, prior to filtering the noise, then add it back in to restore the signal. The algebra 
of this process is presented in Appendix B. We emphasize that W4 is the worst DA for 1/f stripes by at 
least a factor of 3 (Jarosik et al. 2003b), and we limit covariance along scan directions to be <0.1% for all 
other WMAP first-year sky maps. 
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3.3. Pointing and Beam Determination 

3.3.1. Spacecraft Attitude Control and Determination 

The spacecraft attitude is determined from a combination of two autonomous star trackers (ASTs) with 
boresights perpendicular to the spin axis (along the spacecraft ±y axes), two rate gyroscopes, and two digital 
sun sensors. The sensor outputs are combined using a Kalman filter to determine the aspect solution. The 

sensor noise parameters and offsets were initially calibrated in flight during the in-orbit checkout (IOC) 
period in July 2001. By the end of IOC, the final tables were uploaded to the spacecraft. 

Spacecraft quaternions output by the Kalman filter provide the definitive transformation from the 
spacecraft reference frame to the J2000 geocentric inertial (GCI) system. Errors in the attitude solution are 
estimated using the residuals of the individual sensor signals and propagated to the quaternions. After the 
final Kalman filter parameters were loaded, quaternion differences show a noise-like error with a 10" rms. 
In addition to the sensor noise, there is an apparent spin- synchronous error of ~ 10" that is believed to 
be due to propagation errors in the Kalman filter. As discussed below, this error is apparently seen in the 
instrument borcsight determination using Jupiter observations. Since the pointing performance exceeds the 
requirement of 0^9 (root-sum-square for three axes), no correction of the spacecraft quaternions is attempted 
for the first-year processing. Sufficient information exists in the raw telemetry to attempt a correction in the 
future, if it is warranted. Note that random quaternion errors are automatically accounted for in the flight 
beam response maps generated from the Jupiter observations (Page et al. 2003a). 

3.3.2. Instrument Boresight Determination 

As mentioned above, the spacecraft quaternions provide the definitive reference frame for the spacecraft. 
The instrument boresights, 10 each on the A and B sides, are determined from the Jupiter beam maps (Page 
et al. 2003a), which are generated with respect to the spacecraft frame provided by the quaternions. The 
boresight is defined as the location of the peak of a circular Gaussian fit to the main beam. The results of 
this fitting are given in Table 7 as 20 unit vectors in spacecraft coordinates. These are the values used to 
determine instrument pointing in the first-year data processing. The uncertainty in the boresight position 
is ~ 2" per beam in both spacecraft azimuth and elevation. In addition to statistical uncertainty in the 
boresight fits, there are two other potential sources of error in the boresight determination: changes (drifts) 
with time, and errors in the relative time-tagging of quaternion data and science data. 

To test stability, we note that WMAP is in a position to see Jupiter twice per year for about 45 days 
each time. We refer to each ~45 day period as a Jupiter "season". During the first season of each year, 
the boresights scan across Jupiter from roughly ecliptic north to ecliptic south, and vice- versa in the second 
season. As a test of boresight stability, we have generated beam maps from each of the first two seasons 
of data separately, and have fit boresight directions to each. We find the azimuth positions agree to better 
than 3" on both the A and B sides, but the elevation positions differ by 10" on the A side, and a smaller 
amount on the B side. This difference between seasons is consistent with the ^ 10" spin-synchronous error 
in the spacecraft quaternions discussed above. We ignore this small effect in the first year processing, and 
subsume the small systematic error that results into our error estimate for the beam transfer functions, as 
discussed in Page et al. (2003a). 

The relative time-tag accuracy of telemetry packets was tested on the ground. A timing computer 
was set up to simultaneously receive test pulses from both the Attitude Control Electronics (ACE) box 
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and the Digital Electronics Unit (DEU), the two computers that tag the attitude and science data packets, 
respectively. Each of these boxes in turn derives its time from the main "Mongoose" computer on WMAP 
(Bennett et al. 2003a). This test demonstrates a relative time-tag accuracy of 30 /xs between the quaternion 
packets and the science packets. In observing mode, the boresights sweep the sky at a rate of 2?6 s~^, so a 
30 fis time error produces a negligible pointing error of < 0.3" . 

Uncertainty in the spacecraft position is another potential source of boresight determination error. For 
the first year processing we use a predicted ephemcris that is uploaded to the spacecraft approximately 
weekly for on-board use by the Attitude Control System. This solution is returned in telemetry and is the 
basis for the ephemeris data supplied with the first-year release of time-ordered data. The uncertainty in 
these predictions arc < 7 km in position and <1 cm s~^ in velocity (3 cr), relative to the Earth. An error of 
7 km in WMAP's position would result in error of ~2 mas in the apparent position of Jupiter as seen from 
WMAP and is thus completely negligible. 

3.3.3. Beam and Window Function Determination 

Along with gain calibration and noise properties, knowledge of the beam shapes and window functions 
are among the most important aspects of the instrument to characterize for accurate measurements of the 
CMB. Page et al. (2003a) describe in detail the process by which beam maps are formed from in-flight 
observations of Jupiter, and how those maps are transformed to determine the beam window functions. The 
primary result they derive arc a set of 10 beam transfer functions, 6;, one per DA, based on azimuthally- 
averaged beam radial profiles. These transfer functions are included in the first-year data release (Limon et al. 
2003). In addition, they derive a full covariance matrix for each transfer function, Sf;,, which characterizes 
the uncertainty in bi. Typically, the uncertainty for a single DA is about 1-2%, with moderate covariance 
in I. Sec Figure 5 of Page et al. (2003a). These estimates already include a systematic error allowance to 
bound the small range of results obtained from different analysis methods. As described in Hinshaw et al. 
(2003a), the window function covariance matrices are propagated into the Fisher matrix (inverse covariance 
matrix) for the final combined angular power spectrum. Thus the final power spectrum, and the parameter 
fits based on it, already include statistical and systematic window function uncertainties (Spergel et al. 2003; 
Verde et al. 2003; Peiris et al. 2003). 

3.3.4- Effects from Elliptical Beams 

The WMAP beams are moderately elliptical, so the use of azimuthally-averaged radial profiles to de- 
scribe the beam response is an approximation. This approximation is justified in Page et al. (2003a) by 
noting that the WMAP scan pattern produces excellent azimuthal averaging of the beam response in a large 
fraction of the sky. They have placed limits on the variation of the window function across the sky by 
comparing the effective window function in the ecliptic plane, based on a full two-dimensional transform of 
the beam response averaged over the flight range of scan angles, to the fully averaged transform, bi. For 
the three highest frequency cosmology bands, these variations range from 2-3% at Q band to ~1% in V and 
W bands. These variations are consistent with estimates of the angular power spectrum using data at high 
and low ecliptic latitudes separately (Hinshaw et al. 2003a). Since most of the statistical weight at high-Z 
resides in the V and W band data at high ecliptic latitudes, the use of fully averaged beam transforms is 
appropriate, and the systematic error estimate incorporated into 'E^, should encompass any error in this 
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approximation. 

Elliptical beams can also produce errors in the sky maps that are difficult to characterize in a simple 
way. We can define the sky map error due to non-circular beam response as 

where tpbs is the hypothetical noise-free sky map obtained with the actual experimental beam and scan 
pattern, and tdrc is the ideal sky map obtained by convolving the true sky with the averaged beam transform, 

bi. For a diff'erential experiment like WMAP, there are two effects that contribute to Atasym- The first, as 
noted, is incomplete azimuthal coverage in a given pixel, which gives rise to slightly elliptical peak structure 
at low ecliptic latitudes (see below), the second is a localized effect due to echoes from bright Galactic sources 
propagating to other pixels in the map. Specifically, as a bright source is observed in different orientations, 
the differential signal changes with orientation. Since the map-making algorithm must assign one average 
value to the pixel with the bright source, the ring of pair pixels at the beam separation will see an echo with 
a quadrupolar temperature distribution around the ring. We mitigate this effect by incorporating a bright 
source mask in the map- making algorithm, which is invoked as follows. If side A observes a pixel in the 
bright source mask, we only update the sky map accumulator for pixel A, but not for pixel B 

Hobs • tn+l{pA) Hobs • t„+i(pA) + Wi [d{ti) + t„(ps)] nobs(PA) ^ nobs(PA) + Wi 

nobs • t„+l(pB) 7^ nobs • t„+i(pB) - Wi [d{ti) - tn{pA)] nohsiPs) 7^ nobs(Ps) + Wi, (29) 

where the terms are as defined in after equation (19). In this way we obtain an estimate of t{pA), but we do 
not propagate bright echoes to the ring of neighbor pixels, of which pB is one. The mask we use for assigning 
this cut is the same Kp8 processing mask we used for the calibration fits (§2.2.1). 

We have generated a simulation to quantify the errors from both of these effects. Specifically, the 
simulation generates one year of noise-free differential sky signal which includes a model for the flight beam 
ellipticity. We run this data through the flight map- making pipeline to generate sky maps, tobs- We also 

generate convolved maps tdrc using the azimuthally-avcragcd beam transforms appropriate to the beam 
model. The residual map, Atasymj for DA Kl is shown in Figure 13. The K band radiometers have the 
largest beam ellipticity of all the DAs, so this represents a worst case result. The general "mottling" near the 
ecliptic plane results from the relatively limited azimuthal coverage in this region producing elliptical peaks 
and anti-peaks which, in Atasym, are differenced with circular counterparts. This is especially noticeable 
near bright Galactic sources. The rms amplitude of these fluctuations in the Kp2 cut sky is 2 /xK in K 
band, and at least a factor of 2 lower in Q-W bands. The effect of this structure in the power spectrum is 
primarily represented as a variation in the window function across the sky, as discussed above, and in Page 
et al. (2003a) and Hinshaw ot al. (2003b). However, this structure also contributes to the 4-point function 
of the data in the sense that it couples power from different I ranges. This effect is potentially important for 
the interpretation of any gravitational lensing analysis of the WMAP data. 

The flgure also exhibits faint echoes of the brightest sources that evade the map-making cut discussed 
above. We limit localized features in the Kp2 cut sky to less than 10 /uK in K band and less than 2 fiK 
in Q-W bands due to a combination of dimmer sources and more circular beams. We estimate that such 
features occupy <0.1% of the Kp2 cut sky. 
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3.3.5. Far Sidelobe Pick-up 

The WMAP optical system was designed to produce minimal pick-up from signals entering the far 
sidelobes. Barnes et al. (2003) present a complete determination of the WMAP sidelobe response by com- 
bining measurements from a variety of ground-based sources with in-flight measurements of the Moon. They 
produce response maps covering An sr that are included as part of the first-year data release. They then 
use these response maps, with the first-year sky maps, to estimate the systematic artifacts remaining in 
the first-year maps based on the well-justified assumption that sidelobe artifacts are small relative to the 
sky signal. The K band data have the largest sidelobe signal due both to the largest sidelobe spill, and to 
the brightest Galactic signal. The signal was deemed to be large enough, and well enough characterized, 
to warrant a small post-processing correction to the first-year K band map. Limits on remaining sidelobe 
induced artifacts in all the bands are presented in Table 1 of Barnes et al. (2003). 

3.4. Environmental Effects 

3.4- 1- Thermal Effects 

The radiometer gain and offset arc dependent on temperature. There are several aspects of the WAdAP 
design that are critical to mitigating this source of systematic error (Bennett et al. 2003a). The instrument 
is differential, so thermally induced gain changes act on a relatively small offset signal. The observatory 
environment was designed to be as stable as possible, consistent with other goals. For example, all nominal 
thermal control is passive to avoid heaters cycling on and off. The observatory is placed at the second Earth- 
Sun Lagrange point far from the Earth, and the solar panels maintain a fixed 22?5 angle with respect to 
the Sun during normal observing mode. The instrument temperature is monitored with precision platinum 
resistance thermistors (PRTs) to verify the degree to which thermal stability is in fact achieved. 

Temperature variations at the spin period are the most critical since they can induce signals that couple 
relatively efficiently to the sky. However, owing to the relatively fast (129.3 s) spin period and the thermal 
mass of the instrument, any induced signals will have a very red spectrum and thus will couple only to the 
lowest few harmonic modes on the sky (Z<10). In the analysis below we use flight data to estimate the 
susceptibility of the gain and baseline to temperature variations of the instrument. In turn we use limits on 
the instrument's physical temperature variation at the spin period from Jarosik et al. (2003b) to put limits 
on thermally induced artifacts in the time-ordered data, and hence the sky maps. 

Thermally induced signals can enter either through changes in the gain acting on the instrument offset 
or through changes in the offset itself. We show below that the latter are more significant for WMAP. 

The radiometer gain model presented by Jarosik et al. (2003b) describes the gain in terms of the RF 
bias ("total power") of the detector, and the temperature of the FPA. This model tracks thermal variations 
in the gain on the time scale of the RF bias readout (23.04 s), and the map-making algorithm updates the 
gain on this time scale. However, since this sample rate is only a few times per spin, it is possible that a 
systematic temperature variation at the spin period could induce gain changes that are not well sampled by 
this model. As a separate check of gain-induced artifacts, we infer the temperature susceptibility of the gain 
from data taken over a long time period where gain changes are measurable. Results are given in Table 8. 
We combine these measurements with the limits on temperature modulation at the spin period derived by 
Jarosik et al. (2003b) to place limits on gain induced artifacts, as shown in Table 9. We conclude that 
gain-induced signals at the spin period, which might be poorly tracked by the gain model, are <20 nK. 



- 23 - 



The instrument baseline is the product of the gain times the offset. As described in §2.2, we get an initial 
estimate of the baseline from the dipole calibration algorithm. This gives us an estimate of the instrument 
baseline once per hourly precession period. Sample hourly baselines for channels V113 and V114 are shown 
as a function of time over the first year of operation in Figure 14. Also shown is the temperature of the 
instrument FPA over the same time period; there is a clear temperature dependence in the baseline. We 
measure the baseline temperature susceptibility by fitting the hourly baseline estimates to a model of the 
form b(t) = Co + cit + C2ATppA{t) where the Cj are model coefficients and ATppA ~ Tfpa — (Tfpa) is the 
deviation of the FPA temperature from its mean. The most robust susceptibility results come from fitting 
a portion of the data near the time of a partial battery cell failure which occurred on day 2002:054 (GMT) 
(Limon et al. 2003). In response to this event, the spacecraft bus voltage was autonomously commanded 
lower on day 2002:058 (GMT) causing the spacecraft to dissipate less power and thus cool slightly. The 
coefficients are given in Table 8. We have combined the results for the two channels in each radiometer 
because this is the combination that enters into the final sky maps. This has the effect of canceling some 
of the common-mode susceptibility measured in individual channels. As noted above, we combine these 
susceptibility measurements, taken over long time periods, with limits on the temperature variations at the 
spin period (.Jarosik et al. 2003b) to place limits on induced signals at this time scale. The results are given 
in Table 9. We conclude that offset-induced signals at the spin period are <180 nK. 

Slow drifts in the instrument temperature will be largely filtered out by the baseline pre-whitening 
discussed in §2.3.2. The steepest temperature gradient observed during the first year of observation occurred 
just after the above-mentioned battery cell failure. To assess the efficiency of the filtering process, we 
have analyzed the data during this period as follows. We applied the baseline pre-whitening filter to the 
temperature signal, TppAit), to measure how much the cooling gradient was suppressed by the filter. The 
input temperature gradient on day 2002:058 (GMT) was -1.7 mK hr'^ Applying the Kll filter to TppAit) 
yielded an output gradient of —0.1 /xK hr~^, while the W41 filter yielded an upper limit 10 times smaller. 
We conservatively estimate upper limits on residual drift in the ffltered baseline of <10 nK hr~^for the most 
susceptible channels. 

3.4.2. Electrical Effects 

A variable electrical signal on board the observatory could induce an apparent signal in the radiometers. 
Sources of variable electrical signals include the reaction wheels, transponder, bus voltage fluctuations, and 
RF noise coupled to the instrument through the power bus. During the final observatory thermal vacuum 
test, in which the observatory was operating at temperatures close to those achieved in flight, searches for such 
electrically induced radiometric artifacts were conducted (Jarosik et al. 2003a). Upper limits on radiometer 
bus voltage susceptibility, based on ground tests, are given in Table 8. We combine these results with an 
upper limit on bus voltage variations of 3.0 mV rms, measured on-orbit, to conclude that electrically-induced 
signals at the spin period are <40 nK. See Table 9. 

3.5. Miscellaneous Effects 

3.5.1. Radiometer Cross-talk 

A large signal in one radiometer could induce an erroneous signal at the output of another radiometer 
due to electrical cross talk. Such cross talk is not expected but could arise from, e.g., non-ideal amplifier 
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behavior, or other parasitic effects, such as pickup in the wiring harnesses. 

A careful search was made for this effect during the instrument ground tests. Noise diodes were used to 
inject a large signal into one radiometer at a time while the input feeds of all other radiometers were covered 
by absorptive loads. The outputs of the 9 non-driven differencing assemblies were searched for any evidence 
of the injected signal. The tests were run with the amplifiers in the passive channels both on and off in order 
to distinguish pickup mechanisms. No pickup was found in any test. Table 10 gives 2a upper limits on the 
pickup by any DA due to any of the 9 other DAs. The column labeled Electrical gives the results obtained 
from the test with the amplifiers turned off, and the column labeled Radiometric gives weaker limits from 
the test with the HEMTs turned on. The latter limits are weaker because the output noise levels are higher. 

This level of pickup is quite small and could only be of potential concern when WMAP scans across 
Jupiter, the brightest source in the sky for WMAP at L2. The values in the Table 10 are thus referred to 
peak Jupiter signals in each band. For example, the first entry indicates that when Jupiter induces a signal 
of 185 mK in the Wl radiometer, the pickup in the Kl differencing assembly is <30 /iK (95% confidence) 
which is —26.8 dB below the peak Jupiter signal of 14 mK in K band. This signal occurs when the beam of 
the pickup channel is within a few degrees of Jupiter, depending on channel separation in the focal plane, 
and in every case is less than direct radiometric detection of Jupiter in the near side lobes. 

This limit on cross talk implies pickup is completely inconsequential in normal observing mode. Using 
the same Table example, a 100 /xK signal in Wl could cause at most a 200 uK signal in the most susceptible 
the four K band differencing assemblies. 

3.5.2. Source Variability 

Time variable objects are a potential source of contamination for observations of the CMB; see for 
example Sokasian et al. (2001) and references therein. One concern is that an object may grow in brightness 
over the course of WMAP's observations, avoid detection during an initial source survey, and remain un- 
masked during subsequent data analysis. For example, blazars produce relatively rapid and large amplitude 
variability in all wavebands. Long term observations of such objects show that increases in flux by a factor 
of up to a few over a time scale of years can be anticipated (Flett & Henderson 1983; Funis et al. 1982; 
Stevens et al. 1994; Bower et al. 1997). Observations of Zw 2 by Falcke et al. (1999), provide an extreme 
example: a greater than 20-fold increase in brightness, from ^0.1 Jy to ~2 Jy, over a period of less than two 
years. While this object could produce a temperature response of a few hundred /xK in the WMAP data, 
such objects are rare and, if left undetected, would have a minimal effect on cosmological inferences. Tests 
for point source contamination in the WMAP data are given by Bennett et al. (2003c) and Hinshaw et al. 
(2003b). These tests will need to be revisited on an annual basis. 

Another source of concern is that a time variable source in the WMAP data has the effect of broadcasting 
noise to the ring of ^1000 pixels which are separated from the variable source by the beam separation angle 
(6'beam ~ 141°). The point source list derived from the WMAP first-year data is 98% reliable with ~5 
spurious detections at the ^0.5 Jy flux limit of the survey (Bennett et al. 2003c). The nominal point source 
sensitivity of the WMAP telescope is F ^ 200 /xK/Jy, thus a noise level of ~0.1 /uK is expected from variable 
sources that evade detection. 
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4. CONCLUSIONS 

The processing steps used to produce the first-year WMAP sky maps include an initial simultaneous 

estimate of the sky map and the instrument calibration. The instrument gain is then refined using a model 
based on engineering telemetry, and the instrument baseline is refined by the application of a pre-whitcning 
filter. A final archive of calibrated data is produced and used to generate final sky maps using a slightly 
refined iterative algorithm. 

We limit systematic artifacts due to calibration, map-making and environmental disturbances to less 
than ~15 fiK^ in the quadrupolc C2, with tighter limits at higher multipole moments (Table 5). Beam 
transfer functions are measured for each beam with 1-3% over the entire range of multipole moments that 
WMAP is sensitive to (Page et al. 2003a). The covariance matrix of the beam transfer function is propagated 
through to the final power spectrum error matrix. We characterize pixel-pixel covariance matrix and place 
limits on residual stripes in the final maps. 

All major data products from the first year of WMAP observations are being released through NASA's 
new Legacy Archive for Microwave Background Data Analysis (LAMBDA) at http : //lambda. gsfc. nasa. gov/. 

The WMAP mission is made possible by the support of the Office of Space Sciences at NASA Headquar- 
ters and by the hard and capable work of scores of scientists, engineers, technicians, machinists, data analysts, 
budget analysts, managers, administrative staff, and reviewers. We acknowledge use of the HEALPix pack- 
age. 
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A. THE MAPPING FUNCTION 

Equation (4) defines the continuous form of the mapping function, which encodes both the scan strategy 
of an experiment, and convolution due to the beam response. We can relate this to the matrix form, in 
equation (7) as follows. The mapping function evaluated at time U for a finite integration time r may be 
written in terms of the beam response function as 



M{n,U) = - / " dt [aB^(R(i) • n) -0BB(R{t) ■ n)] , (Al) 
where BA{n) is the beam response of the A-side beam, in spacecraft coordinates, normalized to unit integral 

dn^BAin) = l, (A2) 



similarly for the B side, and R{t) is the rotation matrix from sky-fixed (Galactic) coordinates to spacecraft 
coordinates at time t. The terms a and P in equation (Al) account for possible ohmic losses in the A and 
B-side optics that are not necessarily equal (Jarosik et al. 2003b). Since the data are calibrated using the 
modulation of the CMB dipolc, we only need to parameterize the loss imbalance which, following Jarosik 
et al. (2003b), we parameterize as 

O: = 1 -j- Xim, 

/? = (A3) 

Note that loss imbalance is separate from lossless differences in the beam response function, e.g. differences 
in the solid angle of the A and B-side beams. Once the calibration is applied, the differential sky signal is a 
measurement of the form 

At(i) = J dO„t(n) [(1 + Xi^) BA{K{t) • n) - (1 - a;^) BB(R{t) ■ n)] , (A4) 

which still includes the effects of any loss imbalance. We now separately consider how this calibrated 
differential data propagates into the sky maps and the Jupiter beam maps. 

When making sky maps from the calibrated data, each datum is modeled simply as 

At{U) = {1 + Xin,)t{pA) - {I - X-nn) t{pB), (A5) 

where pa is the pixel observed by the A-side beam at time ti, and similarly for p^. That is, each row of the 
mapping matrix in equation (7) has the form 

M{p,ti) = {...,{l+Xin,),...,-il-Xim),...], (A6) 

with non-zero entries in pixel columns pa and ps only. Upon solving for the sky map, this ideally leads to 
an effective beam response of the form 



Bin) = \ 



B^A^ in) , B^\n) 



(1 + Xim) (1 - Xim) 



(A7) 



where is the symmetrized beam response for the A-side beam and similarly for the B-side. For this 
ideal case to obtain, the following must hold: 1) each pixel must be observed equally by the A and B-side 
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beams, which is true to a very good approximation for WMAP, and 2) each pixel must be observed with 
uniform azimuthal coverage. Deviations from these assumptions are considered in the text. 

The beam mapping data is compiled from calibrated observations of the bright source Jupiter. The 
calibrated data taken when side A is observing Jupiter has the form 



At 



j rfOn t j(n) (1 + Xi^) BaCR ■ n) + At^ky , (A8) 



where tj(n) is the brightness temperature of Jupiter in the direction n and Atsky is the background sky 
temperature difference, which is subtracted during processing. An analogous equation holds for the B-side 
data. Assuming the beam response is constant over the extent of Jupiter, the integral reduces to 

At = TjQj (1 + Xim) BaCR ■ nj), +Atsky, (A9) 

where Tj is the disk brightness temperature of Jupiter, and flj is its solid angle. Beam maps are compiled 
by binning the corrected data At — Atgky as a function of n, in spacecraft coordinates. This produces maps 
proportional to the beam response 

TA{n) = TjQj{l + Xi^) BA{n), (AlO) 

and similarly for the B-side. 

Ultimately, we wish to compute the transfer function of the symmetrized beam response. This may be 
obtained from the symmetrized beam maps as 



B{n) ^ 



(l+a;im) (1 - a;im) 



(All) 



We don't know the brightness temperature of Jupiter a priori, but since this is an overall normalization 
factor, we are free to normalize the final transfer function to 1 at ^ = 0. 



B. MAP-MAKING WITH FILTERED DATA 

§2.3.2 presents the filtering algorithm used to determine the final instrument baseline. This process 
included an estimated sky signal subtraction based on the initial sky maps produced with the hourly cal- 
ibration. In this Appendix we derive the noise properties of sky maps produced with this filtered data. 
In the following, we assume the time-ordered data has a noise covariance N = (nn^) that includes a 1// 
component, and that we have a pre- whitening filter F such that 

((Fn)(Fn)'^) ocl. (Bl) 



B.l. Map- making with filtered signal -|- noise 

We could filter the full data prior to any calibration or sky map estimation, then deconvolve the effects 
of the filter in the subsequent data processing. The input data would have the form 

d' = Fd = FMt-FFn. (B2) 
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Then, in order to obtain an unbiased sky map estimate, we would need to evaluate the sky map estimator 

t' = (M'^F'^FM)-i • (M'^F^d'), (B3) 

which deconvolves the action of the filter on the sky signal. Since the 1// noise in the WMAP data is 
relatively small, implementing this estimator for the first-year sky maps was deemed unnecessary and would 
likely have delayed the release of the maps. The alternative is to filter only the noise by subtracting an 
estimate of the sky signal prior to filtering, then adding it back in to the time-ordered data prior to making 
new maps. 



B.2. Map-making with filtered noise 

Let to be the sky map estimated from unfiltered data, using the hourly calibration. This is related to 
the true sky signal by 

to = Wd (B4) 
= W-(Mt-hn) (B5) 
= t + Wn (B6) 

where W = (M^M)~^ • is the map-making operator defined in §2.1, and we have used the fact that 
W • M = 1. We use this sky map to subtract a sky signal from the time-ordered data prior to filtering, then 
we add it back in after filtering. This produces a filtered data set 

di = F • (d - Mto) + Mto (B7) 
= F • (Mt + n - Mt - MWn) -h Mt -h MWn (B8) 
= Mt -h Fn + (1 - F) • MWn. (B9) 

This time scries data set consists of an unbiased sky signal Mt, a white noise term Fn. and a residual noise 
term (1 — F) • MWn which is due to the off-diagonal "wings" of the filter (1 — F) acting on the noise from 
the initial sky map estimate, MWn. 

We can make a map ti from the data di using the algorithm of §2.1 

tl = Wdi (BIO) 

= t + WFn + W- (1 -F) -MWn (Bll) 

= t + WFn + (1 - WFM) • Wn (B12) 

= t + WFn + R-Wn, (B13) 

where we have again used W • M = 1 and wo have defined a "residual" operator R = (1 — WFM) which is 
small in the sense that only off-diagonal terms in F contribute to it. This is most easily seen if we recall that 
the filter operator is 1 on the diagonal and has small off-diagonal terms. We can then write F = 1 — E from 
which it follows R = WEM. It follows that ti is an unbiased estimate of t that includes a white noise term 
WFn (this noise still contains the small beam separation covariance), and a residual noise term, R • Wn, 
due to the noisy sky signal estimator used in the filtering process. This latter term is the "excess" noise seen 
in the W band single DA maps after filtering (Hinshaw et al. 2003a). 
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The residual noise term can be reduced somewhat by iterating the filter algorithm a second time, using 
ti as a sky signal estimator. It is straightforward to show that, after some algebra, the resulting map is 

t2 =t + WFn + R-WFn + C»(R2), (B14) 

where we have neglected a term of order • Wn. The residual noise is reduced slightly since R is now 
acting on the white noise WFn instead of the full noise Wn. But since the white noise dominates, this 
is a relatively insignificant improvement. It is also clear that subsequent iterations of the filtering only 
contribute higher order corrections which are negligible. This convergence has been verified with the flight 
data. Another feature wc have verified with the flight data is the fact that the excess noise term decreases 
with time. The reason for this is simply that R = (1 — WFM) = WEM gets smaller with additional data 
because the map-making operator W gets smaller as more observations accumulate in its "denominator", 
(M'^'M)"^. The first-year maps were only filtered once, as per equation (B13), because the improvement in 
noise properties was not deemed sufficient to warrant the additional processing time. 



C. MAP-MAKING WITH POLARIZATION 

WMAP observes the sky with two orthogonal linear polarization modes per feed, thus it is sensitive to 
the 3 Stokes parameters I, Q, and U. This Appendix outlines the algorithm with which these parameters can 
be mapped using the differential data from WMAP. The approach is an extension of the iterative method 
in §2.1.1 introduced by Wright et al. (1996). 



C.l. Polctrization Mapping with Total Power Data 

Suppose we observed the sky with a single beam, total power radiometer that is sensitive to a single 
linear polarization, denoted mode ^1. In terms of the Stokes parameters, the temperature observed by the 
instrument at time t in pixel p would be 

di (i) = i(p) + q(p) cos 27 + u(p) sin 27, (CI) 

where i, q, and u arc Stokes parameter maps in units of temperature, and 7 is the angle between the 
polarization axis of the beam and the chosen reference direction for pixel p (the choice of reference direction 
is discussed below). Note that we adopt the common convention that the instrument response reduces to 
the total temperature in the limit of an unpolarized source. 

The signal in the orthogonal polarization channel, which is fed by the other port of the ortho-mode 
transducer (GMT) and denoted mode #2, is given by 



d2{t) = i(p) + q(p)cos2(7+-) + u(p)sin2(7+-) (C2) 

= i(p) — q(p) cos 27 — u{p) sin 27. (C3) 
We can isolate the intensity and polarization signals by taking sums and differences 

d{t) ^ i(di + d2) = i(p) (C4) 

p{t) = ^(di-d2) = q(p)cos27 + u(p)sin27. (C5) 
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Given noisy data, we can estimate the intensity i{p) by averaging all the data d{ti). For the polarization we 
can only estimate a linear combination of q and u from a single observation. However, if we have several 
observations of pixel p with a variety of polarization angles 7, we can estimate q and u in a given pixel by 
minimizing x^, defined as 

2 _ [P(*i) - ^{P) COS 27i - u{p) sin 2'yif 

X =2^ , i»^oj 

where i is a time-ordered data index and the sum is over observations within pixel p, and 7j is the polarization 
angle for the i*'' observation. The best-fit values for q and u are given by 



(C7) 



where Ci = cos27i, si = sin27i, and A = Ei "^f /'^f Ei /'''f ^ (Ei '^j*'/'''?)^ the determinant of the 
normal equations matrix. The standard errors for q and u are given by the inverse of the normal equations 
matrix 

< = ^E^i/'^l (C8) 

i 

= (C9) 

i 

< = (CIO) 

In the limit of uniform azimuthal coverage and constant noise per observation (cr, = ctq), the matrix elements 
in the linear system reduce to 



EiCf E^CiSi \ f N/2 



HiCiS. EiSi J V N/2 



(cii) 



where A'^ is the number of observations of pixel p. In this limit, the noise in q and u is equal and uncorrelated 
and reduces to 

O-q = CTu > d — (Jo- (C12) 

Thus the noise in each polarization component is times noisier than in the intensity i. 



C.2. Polarization Mapping with Differential Data 

We now generalize to the case of polarization mapping with differential input data. For clarity, we 
first consider the case where the loss in the two sides A and B are equal. We generalize to the case with 
unbalanced loss in the next subsection. In the case of WMAP, one radiometer in a differencing assembly 
(radiometer #1) is fed from the axial port of the OMT and the other (#2) is fed by the lateral port (Jarosik 
et al. 2003a). Following equations (C1,C3) the differential signals from radiometers 1 and 2 are 



di{t) = -^{di3 + du) = i{pA) +<l{PA)cos2^A + u{pA)sm2^A (C13) 
- {{pb) - q(pB) cos 27b - u{pb) sin 273, (C14) 
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and 

^2(0 = ^(d23 + d24) = Kpa) - q(PA) COS - u{pa) sin 27a (C15) 
- i{pB) + (i{pb) cos 2 jB + u{pB)sm2jB, (C16) 

where 7a is the angle between the axial polarization plane and the reference direction in the pixel seen by 
the A beam, and similarly for 7^. We take sums and differences of the two signals to isolate the unpolarized 
and polarized portions of the signal 

d{t) = l{di+ da) = i{pA) - i{pB) (C17) 

p{t) = ^(di - da) = q(pA)cos27A +u(pA)sin27A - q(pB)cos27s - u{pB)sm2jB- (C18) 

An iterative scheme for making maps of q and u follows the form used for intensity maps (§2.1.1). 
Let q,i and u„ be the n^^ estimates of q and u respectively. Estimates for qn+i and u„+i are obtained 
by combining the per-pixel fitting algorithm in equation (C6) with the iterative map-making algorithm, as 
follows 

^2 ^ ^ [P'jti) - qn+i(p) cos27, - Un+i{p) sin27i]^ ^^^^^ 

where the sum is over all observations of pixel p by either the A- or B-side beam, and p'{t) is the polarization 
data corrected with an estimate of the signal in the opposite beam 

/.^ ^ f P{ti) + qnbs) cos 27b + u„(pb) sin 275 beam Aep 
~ I -p{ti) + cin{pA) COS 27a + u„(pa) sin 27a beam B e p. 

The best-fit solution for q„+i and u„+i is then 



U„+i(p) J A V - Y.^C^sJaf Y.^cf/af J \ Y.iSiP'iU)/<7^ 



(C21) 



where the sum on i is as defined above. The uncertainties are as given for the total power case, where ai is 
now the uncertainty per differential observation, p{ti). 



C.3. Map Making w^ith Unbalanced Differential Data 

We now generalize to the case of map-making with unbalanced differential input data. In this case, 
losses in the A and B-side front ends are different and the differential signal is of the form 

di(i) = (l-Fa;i„,,i)[i(pA)+s(pA,7A)]-(l-a;ix„,i)[i(ps)-Fs(pB,7B)] (C22) 
d2(i) = {1 + Xiin,2) IKpa) - s(pA, 7^)] - (1 - a;iin,2) [i(ps) - s(pb, 7b)] > (C23) 

where aJim.i, a;im,2 are the loss imbalance parameters in radiometers 1 and 2, respectively, as defined in Jarosik 
et al. (2003b), and s(p, 7) is short-hand for the linear combination of Stokes parameters 



s(p, 7) = q(p) cos 27 -|- u(p) sin 27. 



(C24) 
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As before, we take sums and differences of the two signals to isolate the unpolarizcd and polarized portions 
of the signal. First, define the mean imbalance and the "imbalance in the imbalance" as 

Sim = ^(a;im,i + a;im,2) (C25) 
Sxiyn = ^(a;im,i - a;im,2), (C26) 

then 

d{t) = ^(di + d2) = i{pA) - Kpb) + Sim [Kpa) + Kpb)] + Sxi^n [s(pa,7a) + s(ps,7b)] 

= (1 + Sim) Kpa) - (1 - Sim) Kpb) + 0{6Xi^) (C27) 
p{t) = ^(di - d2) = s(pa,7a) - s(pb,7s) + Sim [s(pa,7b) + s(pb,7s)] + ^a;im [Kpa) + Kpb)] 

= (1 + Sim) S(PA, 7a) - (1 - Sim) s(pb, 7b) + C(feim)- (C28) 

Note that the term of 0{5xini) in equation (C27) is negligible, because i s, but the term in equation (C28) 
must be considered more carefully. First note that Sxi„i is small - from Table 3 in Jarosik et al. (2003b), the 
largest value is 0.35% in W2, with most values being ^0.1%. Second, the multiplier, (i(p/i) + Kpb)) does 
not modulate with polarization angle, 7, thus it is effectively an offset term that is further suppressed by 
the map-making demodulation. Finally, the term is out of phase with the differential signal (i(pA) — Kpb}) 
so it does not effectively couple to the sky. The effect of this term in the first-year data is further considered 
by Kogut et al. (2003). 

We generalize the differential map-making algorithm to account for loss imbalance as follows. For 
intensity, let i„ be the n^^ estimate of i, then 

i.PB)]/il + ) beam A G p //-ion\ 
° - 1 [-d{U) + (1 + Sim) Upa)] /(I - Sim) beam B G p. 

The updated intensity map is then estimated by binning the corrected data 

= E^^ld^M, (C30) 

where Wi is the normalized statistical weight of each observation 



(1 + .Tim)^(To/'^j^ beam A Gp 
(1 - Sim)^crg/fT,^ beam B e p. 



For polarization, let q„ and u„ be the n*^ estimates of q and u respectively, and define s„(p, 7) = 
q„(p) cos 27 + u„(p) sin 27. Estimates for qn+i and u„+i are obtained by combining the per-pixel q and u 
demodulation with the iterative map-making algorithm. We define as follows 

2 _ [P'(*i) - Qn+iW cos27i - u„+i(p) sin27j]^ 

where the sum is over all observations of pixel p by either the A- or B-side beam, and 

'u.) = l [+p(^i) + (l-Sim)sn(pB,7B)]/(l + Sim) beam A Gp 
~ \ [-p{ti) + (1 + Sim)s„(pA, 7^)] /(I - Sim) beam B G p. 
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The best-fit solution for q„+i and u„+i is then 

<ln+l{p) \ _ ^ f Y^iWisf -Y^iWiCiSi \ ( YiWiCip'iU) 



where the weight Wi is the same as equation (C31). 



(C34) 



C.4. The Choice of Reference Direction 

Wc adopt the sign convention of Zaldarriaga & Seljak (1997) in whieh the polarization components are 
defined in a right-handed coordinate system with the z-axis pointed outward towards the slcy. The Stol?;es 
parameters are defined with respect to a fiducial direction in each pixel on the sky. We adopt the convention 
in which the reference direction is aligned with the local Galactic meridian, i.e., the great circle connecting 
a given point to the North Galactic Pole. The unit vector, n, tangent to this great circle, pointing in the 
direction of the North Pole, is given by 

n = I X e = I X — 

sm 

where I is the unit vector along the line of sight of the current observation, e is a unit vector pointing east 
from I, z is the unit vector to the North Galactic Pole, and is the polar angle (co-latitude) of I. 

For reference, we give formulae for computing the factors cos 27 and sin 27 here. Let I be the unit vector 
along the line of sight, w be the unit vector pointing west from I, nhe the unit vector pointing north from 
I (the polarization reference direction) and p be the unit vector along the polarization plane defined by the 
axial port of the OMT. Then, for both the A and B sides, we have 

cos 7 = p - h 

sin 7 = p- u) 
cos 27 = 2 cos^ 7 — 1 
sin 27 = 2 sin 7 cos 7. 

Note that this defines a right-handed coordinate system with {x, y, z) axes (n, w, I) whose z axis is 
oriented outward following the sign conventions of Zaldarriaga & Seljak (1997). 
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Fig. 1. — A residual sky map, tout — tin, from an "ideal" one- year simulation of Q2 data, designed to test 
the iterative map-making algorithm presented in §2.1. The input sky map included realistic CMB signal 
with a peak-to-peak amplitude of ~ ±420 /iK, and a Galactic signal with a peak brightness of ^50 mK. 
The rms structure in this map is <0.2 /iK, after accounting for the 0.15 /xK noise that was introduced to 
the simulation to dither the digitized signal. The map is projected in ecliptic coordinates and shows the 
anisotropy mode that is least well measured by WMAP, due to a combination of the scan pattern and the 
beam separation angle. This residual level is the result of 50 iterations of the algorithm - more iterations 
would reduce it even further. 
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Fig. 2. — Convergence of the dipole-based gain solution for a selected WMAP radiometer channel (K113) 
based on a one-year simulation. This simulation was generated with an input gain of 300.0 du mK~^, and 
minimal noise. The first iteration, which assumes the sky model has only a dipole component, leaves residual 
gain errors of up to 7%, due to the projection of the relatively bright Galactic emission onto the dipole model. 
After 30 iterations of the simultaneous fit described in §2.2.1, the residual errors in the gain solution arc less 
than 0.1%. 
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Fig. 3. — The hourly gain and basehne fit described in §2.2.1 from the flight data for channels K113 and 
V113. The top two panels show the gain solution, the bottom two the baseline. Note that the gain is stable 
to ~5% over the first year (see also Table 3). The variable noise is due to the changing projection of the 
scan pattern on the CMB dipole over the course of a year. The instrument baselines have a typical drift of 
5-10 mK over the year. The channel V113 exhibits one of the clearest thermal susceptibilities of the WMAP 
radiometers, though we show in §3.4.1 that the induced systematic signal is negligible. See also Figure 14. 
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Fig. 4. — The left-hand panels show the measured auto-correlation function, C{At)/C{0), for selected 
radiometers, of WMAP time-ordered data, after subtracting a model sky signal based on the; initial sky 
maps. The model fits are indicated by an x at a lag of 1 observation and by straight lines for At > t 
(see §2.3.2). AU of the WMAP DAs except W4 have a covariance of <1% at non-zero lag (Table 4). The 
anti-correlation at lag -^2000 s is due to the subtraction of the hourly baseline as a pre-filter. The middle 
panels show the prc-whitcning filter, in the frequency domain, that is applied to the time-ordered data after 
a model sky signal has been subtracted. The vertical dashed line indicates the spin frequency, 7.7 mHz, and 
the number indicates the fraction of power transmitted by the filter at the spin frequency. The right-hand 
panels show the measured auto-correlation function for selected channels of WMAP time-ordered data, after 
pre-whitcning, on the same scale as the left panels. The apparent change in noise level at different lags in 
C{At) is due to a step- wise change of bin size in At. 
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Fig. 5. — The time-ordered data for channel W424 before (black) and after (red) applying the pre- whitening 
filter. The data are boxcar averaged over a 46.08 s window to show the low frequency noise in the unfiltered 
data. Without averaging, the data before and after filtering are virtually indistinguishable. Note that baseline 
variations in this channel are of order 2 mK on a one-hour time scale, as expected given the measured 1// 
knee frequency of this radiometer (Jarosik et al. 2003b). W4 is the worst differencing assembly from the 
standpoint of 1// noise. 
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Fig. 6. — Simulated and flight difference maps showing the structure that is removed from the maps by the 
pre-whitening filter. All four maps are differences between sky maps generated before and after baseline 
filtering. The maps are projected in Galactic coordinates and the temperature scale on each is ±50 /iK. The 
"blobs" of white noise along the ecliptic plane can be ignored. They arise from differences in the handling 
of planet flags in the two forms of the map-making code. The top two panels show W3 and W4 data from 
a one-year simulation that includes flight- like 1// noise in the time-ordered data. The bottom two panels 
are the same for the flight W3 and W4 maps. Note the very different structure between W3 and W4, due 
to different 1// knee frequencies (Jarosik et al. 2003b). Note also that the simulation captures the basic 
structure of the flight data very well. 
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Fig. 7. — The recovered gain solutions for channel VI 13 in a flight-like simulation (top) and in the flight 
data (bottom) . The "noisy" black traces show the hourly baseline binned in 24-hr samples (to reduce noise) 
and the green traces are the best fit gain model (§2.3.1). For the simulation, the input gain used to generate 
the simulated data is shown in grey. In the simulation, the absolute gain is recovered to better than 0.1% 
in all 40 channels, and the binned hourly gain is everywhere within ~0.2% of the gain model, and the input 
gain. Gain changes are well tracked by the pipeline. 




Fig. 8. — Angular power spectra of the 10 resdidual maps tout — tin generated from the flight-like one-year 
simulation. In each case the spectra were evaluated in the Kp2 cut sky (Bennett et al. 2003c). Table 5 
quantifies structure in these maps beyond flat white noise. Note that most features are restricted to 1^25 
but with an amplitude that is still much less than the sky signal in this range. The residual effects of I// 
noise are seen in the gradual rise of the noise spectrum at low / in W4. See Hinshaw et al. (2003b) for further 
discussion of this. 




Fig. 9. — Sum and difference maps generated from the flight Q, V, and W band data, as indicated. To 
reduce the noise, all maps have been binned in larger pixels (HEALPix Nside — 64) and displayed with a 
temperature scale of ±100 /^K. As discussed in the §3.2, the only apparent structure in the difference maps 
is due to residual Galactic contamination owing to the fact that the effective frequencies of the DAs are 
slightly different. This does not affect signals with a CMB spectrum, because the calibration source (the 
CMB dipole) has the same spectrum. See Bennett et al. (2003b) for higher resolution images of the signal 
maps. 




Fig. 10. — Two-point correlation functions of A12 difference maps for three different DAs. With the exception 
of a ~0.3% bUp at the beam separation angle, 6'beam ~ 141°, there is relatively little structure in the difference 
maps (see §3.2). The two-point functions of these maps provide a good representation of the angle-averaged 
pixel-pixel noise covariance in the flight maps. 
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Fig. 11. — Angular power spectra of signal and noise maps for each DA. In each panel, the upper red and 
green traces are the spectra of the null maps, A34 and A1234, respectively. The lower green trace is the A12 
map, and the black trace is the signal map. The blue curve is our best estimate of the underlying CMB 
signal from Hinshaw et al. (2003b). The pairing of white noise levels is discussed in §3.2, Table 6 presents a 
measure of structure in the difference spectra, which are remarkably flat. 
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Fig. 12. — An estimate of pixel-pixel noise covariance in a W3 and W4 noise map along the scan direction 
before and after filtering (top and bottom, respectively). See §3.2 for a description of the processing steps 
used to produce these data. The stripe covariance is negligible in W3, and <0.2% in W4 for lags >0.1 s 
(pixel separation >0?25). All other DAs will have at least 2-3 times lower covariance than W4. 



Fig. 13. — Residual map from a K band elliptical beam simulation. The output map was generated from 
a one-year simulation of data with an elliptical beam response. The residual map shown was generated 
by subtracting the underlying sky signal convolved with the nearest effective circular beam response. This 
remaining structure contributes to the four-point fluctuation spectrum. The scale of the color range is ±10 
fiK. The rms structure in the Kp2 cut sky is 2 fiK. See §3.3.4. 
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Fig. 14. — An illustration of the mild thermal susceptibility of the instrument baseline. The top panel 
shows the temperature of the instrument Focal Plane Assembly (FPA) over the course of the first year. The 
second and third panels show the hourly baseline solution for channels V113 and V114, which are among the 
most thermally susceptible. Note that the thermal baseline response is mostly common-mode. The channel 
combination that contains sky signal is the difference between channels 3 and 4, thus most of this response 
cancels. On day 2002:054 (GMT) a partial battery cell failure led to a commanded decrease in spacecraft 
bus voltage with a corresponding decrease in overall power dissipation and spacecraft temperature. This 
event provides a clean measurement of the instrument baseline thermal susceptibility - the bottom panel 
shows a close-up of the V114 baseline near this event. The dashed line is a fit to a model including a term 
proportional to dh/dTppA. The best-fit susceptibility values for all channels are given in Table 8. See Limon 
et al. (2003) for a complete discussion of WMAP's first-year thermal profile. 



Table 1. WMAP data processing notation 



Synibol(s) Description 



t, ti Time, time of iih observation, in s 

r Integration time per observation, in s 

Np Number of pixels in a map, — (Np — 1) 

Nf Number of time-ordered data points 

p HEALPix pixel number 

Pa, Pb a & B-sidc pixels at time t 

Oheam Separation of A & B-side beams, in degrees 

7 Polarization angle 

7a, 7b a & B-side polarization angles, at time t 

Ci, Si, CiSi cos 27, sin 27, cos 27 • sin 27, at time ti 

t{p) Sky map, in mK 

i{p) Sky map. Stokes parameter I, in mK 

q(p), u{p) Sky map, Stokes parameters Q, U, in mK 

t(p), t„(p) Estimated sky maps, in mK 

nobs(p) Number of observation of pixel p 

tc{p) CMB anisotropy map, in mK 

tg (p) Galactic foreground map, in mK 

ts{p,t) Time-dependent source map (Sun, Earth, Moon), in mK 

At(t) Time-ordered differential sky signal, in mK 

Atd(t) Time-ordered COBE dipole signal, in mK 

Atv(t) Time-ordered local velocity dipole signal, in mK 

Ata(t) Time-ordered anisotropy signal, At — Atd, in mK 

c{t) Time-ordered raw data, single channel, in du 

dij{t) Time-ordered data, radiometer i, channel j, in mK 

d{t) Time-ordered intensity data, co-added channels, in mK 

p{t) Time-ordered polarization data, co- added channels, in mK 

g(t) Rad. responsivity (oc gain), single channel, in du mK~^ 

n{t) Rad. noise, single channel in calibration, 

4 co-added channels in map-making, in mK 

h{t) Rad. baseline, single channel, in du 

Qk, bk Hourly gain, baseline, single channel, k*^precession 

CTj rms noise, ^''^observation, single channel in calibration, 

4 co-added channels in map-making, in mK 

CTo Mean rms noise, single or co-added channels, in mK 

a;im, Ximj Loss imbalance parameter (radiometer j) 

C{At) Auto-correlation function of noise, in mK^ 

C\, A, B Auto-correlation function model parameters 

w{f) Pre-whitening filter, Fourier space 
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Table 1 — Continued 



Symbol(s) Description 



F Pre-whitening filter, time domain, Nt x Nt matrix 

N Tirnci-ordcrcid noise covariancc, Nt x Nt matrix, in mK^ 

M Mapping function, Nt x Np matrix 

W Map-making operator, (M'^M)-^ • M^, Np x Nt matrix 

S Pixel-pixel noise covariancc. Np x Np matrix, in mK^ 

D Reduced inverse noise, (M^M) = ctq S^^, Np x Np matrix 

Ajj Difference map from channel combination ij 
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Table 2. WMAP Attitude Control System Requirements 



Parameter 


Requirement 


Performance 


Precession rate (d^/dt) 


-0?ls-i ±6.3% 


-0?ls-i ±3.6% 


Spin rate {dijj/dt) 


2? 784s- 1 ± 5% 


2?784s-i±0.13% 


Sun-spin angle (9) 


22?5 ± 0?25 


22?5±0?023 



Table 3. WMAP Dipole-Based Gain Summary 



Channel 3 Channel 4 



Radiometer 


{9k) 






{9k) 








(du mK-i) 


(%) 


(%) 


(du mK-i) 


(%) 


(%) 


Kll 


-0.974 


7.4 


0.66 


+0.997 


6.8 


0.66 


K12 


+ 1.177 


6.2 


0.75 


-1.122 


6.4 


0.75 


Kail 


+0.849 


4.7 


0.75 


-0.858 


5.0 


0.75 


Kal2 


-1.071 


5.1 


0.75 


+0.985 


5.3 


0.75 


Qll 


+1.015 


4.5 


0.94 


-0.948 


4.4 


0.94 


Q12 


+0.475 


5.3 


1.03 


-0.518 


5.5 


1.03 


Q21 


-0.958 


5.8 


0.94 


+0.986 


6.0 


0.94 


Q22 


-0.783 


2.9 


1.22 


+0.760 


2.8 


1.22 


Vll 


+0.449 


4.5 


1.50 


-0.494 


4.5 


1.50 


V12 


-0.532 


4.0 


1.40 


+0.532 


4.7 


1.40 


V21 


-0.450 


4.7 


1.22 


+0.443 


5.1 


1.22 


V22 


+0.373 


3.2 


1.59 


-0.346 


3.0 


1.59 


Wll 


+0.311 


5.1 


2.25 


-0.332 


4.1 


2.43 


W12 


+0.262 


3.5 


2.62 


-0.239 


6.0 


2.71 


W21 


-0.288 


4.6 


3.09 


+0.297 


3.8 


2.53 


W22 


+0.293 


6.1 


2.43 


-0.293 


6.3 


2.62 


W31 


-0.260 


3.3 


2.25 


+0.281 


3.8 


2.34 


W32 


-0.263 


3.6 


2.62 


+0.258 


3.4 


2.43 


W41 


+0.226 


6.0 


4.40 


-0.232 


5.7 


4.21 


W42 


+0.302 


6.3 


3.28 


-0.286 


5.9 


3.37 



''Peak-peak variation in the daily mean gain, indicates range of gain drift 
during first year. 

''Mean statistical uncertainty per hour. 



- 52 - 



Table 4. Auto-correlation Model Parameters'' 



Radiometer 






A 




B 


"'(/spin)'' 




— U.UUOo 


n 

u. 


nni 1 
.uu± ± 


n 

u, 




u.yoo 


K12 


0.0008 


0, 


.0011 


0, 


.00040 


0.963 


Kali 


-0.0075 


0, 


.0015 


0, 


.00048 


0.960 


Kal2 


-0.0031 


0, 


.0006 


0, 


.00019 


0.984 


Qll 


0.0044 


0, 


.0018 


0, 


.00053 


0.934 


Q12 


-0.0088 


0, 


.0007 


0, 


.00024 


0.978 


Q21 


0.0124 


0, 


.0088 


0, 


.00282 


0.754 


Q22 


0.0178 


0, 


.0128 


0, 


.00415 


0.686 


Vll 


0.0010 


0, 


.0001 


0, 


.00005 


0.989 


V12 


0.0034 


0, 


.0014 


0, 


.00048 


0.925 


V21 


-0.0038 


0, 


.0010 


0, 


.00033 


0.951 


V22 


0.0087 


0, 


.0093 


0, 


.00320 


0.689 


Wll 


0.0158 


0, 


.0062 


0, 


.00211 


0.680 


W12 


0.0048 


0, 


.0005 


0, 


.00019 


0.950 


W21 


0.0207 


0, 


.0071 


0, 


,00262 


0.644 


W22 


0.0167 


0. 


.0053 


0. 


.00187 


0.701 


W31 


0.0062 


0, 


.0006 


0, 


.00021 


0.943 


W32 


0.0077 


0, 


.0002 


0, 


.00007 


0.975 


W41 


0.0562 


0. 


.0323 


0, 


.01152 


0.374 


W42 


0.0393 


0, 


.0194 


0, 


.00692 


0.461 



''See equation 26 for model definition. All parameters 
are dimensionless. 

^Derived filter response at the spin frequency, 7.7 mHz. 
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Table 5. Calibration and Map-Making Error Limits'^ 



DA 


C2 

(mK^) 


(C;)3_io 
(mK2) 


(C';)ii_ioo 


(mK2) 




(mK^) 


Kl 


-21.4 


0.6 


0.08 


42.9 


1.1 


0.03 


Kal 


18.5 


1.3 


0.06 


37.0 


2.5 


0.01 


Ql 


59.6 


1.2 


0.14 


118.9 


2.2 


0.01 


Q2 


7.3 


0.9 


0.13 


14.4 


1.6 


0.02 


VI 


3.9 


0.6 


0.21 


7.4 


0.7 


0.01 


V2 


-6.1 


0.8 


0.19 


12.6 


1.2 


0.03 


Wl 


-2.6 


1.4 


0.49 


6.0 


2.0 


0.10 


W2 


12.0 


0.7 


0.62 


22.9 


0.4 


0.15 


W3 


4.3 


0.4 


0.65 


7.3 


0.4 


0.07 


W4 


-6.6 


3.3 


0.90 


14.5 


5.4 


0.55 



''AH values derived from a one-year simulation of WMAP data. The first 
3 data columns give the mean power in the residual map tout — tin from the 
simulation. The last 3 columns give an estimate of the systematic error due to 
calibration and map-making, as defined in §3.1. For comparison, the average 
power in the CMB in each band is C2 ^ 130 /xK^, (Cj)g_-,^Q ~ 150 ^K^, and 
(t^;)ii-ioo~6MK^ 
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Table 6. Difference Map Statistics 



DA 


Diff. map* 


Cf6'hp«Tn)/CfO) 

^ V Deam / / ^ V / 


Co 


\^'/3-10 


\^'/ll-100 


I2 


1 I3~10 
(wK2) 

(^pj-v j 


1 Ill-lOO 

(wK2) 

(^piv j 


Kl 


Ai2 


0.160 


107.30 


1.77 


0.12 


107.25 


1.72 


0.071 


Kl 


A34 


0.014 


6.43 


0.11 


0.06 


6.37 


0.05 


0.003 


Kl 


A1234 


0.030 


13.56 


0.30 


0.07 


13.50 


0.24 


0.011 


Kal 


A12 


0.0057 


2.11 


0.14 


0.06 


2.06 


0.08 


0.002 


Kal 


A34 


0.0022 


0.11 


0.08 


0.07 


0.05 


0.01 


0.002 


Kal 


A1234 


0.0028 


0.01 


0.09 


0.07 


0.05 


0.03 


0.003 


Ql 


A12 


0.0035 


1.08 


0.14 


0.10 


0.98 


0.04 


0.004 


Ql 


A34 


0.0032 


0.49 


0.16 


0.13 


0.37 


0.03 


0.003 


Ql 


A 1234 


0.0044 


0.31 


0.17 


0.13 


0.19 


0.05 


0.003 


Q2 


A12 


0.0031 


0.09 


0.14 


0.09 


0.00 


0.05 


0.005 


Q2 


A34 


0.0030 


0.28 


0.13 


0.10 


0.18 


0.03 


0.001 


Q2 


A1234 


0.0025 


0.03 


0.11 


0.10 


0.07 


0.01 


0.004 


VI 


A12 


0.0038 


4.40 


0.35 


0.16 


4.24 


0.19 


0.006 


VI 


A34 


0.0032 


0.22 


0.35 


0.28 


0.05 


0.09 


0.016 


VI 


A1234 


0.0024 


0.26 


0.30 


0.27 


0.01 


0.04 


0.006 


V2 


A12 


0.0043 


1.72 


0.15 


0.13 


1.59 


0.02 


0.006 


V2 


A34 


0.0026 


0.57 


0.25 


0.22 


0.35 


0.03 


0.003 


V2 


A 1234 


0.0033 


1.05 


0.38 


0.23 


0.83 


0.16 


0.008 


Wl 


A12 


0.0036 


5.45 


0.36 


0.36 


5.11 


0.02 


0.025 


Wl 


A34 


0.0033 


2.10 


0.74 


0.56 


1.56 


0.21 


0.019 


Wl 


A 1234 


0.0034 


0.14 


0.75 


0.57 


0.39 


0.22 


0.037 


W2 


A12 


0.0029 


1.20 


0.51 


0.44 


0.79 


0.10 


0.021 


W2 


A34 


0.0026 


0.12 


0.57 


0.57 


0.44 


0.00 


0.009 


W2 


A 1234 


0.0028 


0.26 


0.56 


0.59 


0.31 


0.01 


0.025 


W3 


A12 


0.0035 


4.95 


0.49 


0.47 


4.49 


0.03 


0.013 


W3 


A34 


0.0031 


8.84 


0.90 


0.75 


8.12 


0.19 


0.039 


W3 


A1234 


0.0039 


2.91 


1.11 


0.72 


2.20 


0.40 


0.007 


W4 


A12 


0.0030 


1.50 


0.44 


0.48 


1.05 


0.00 


0.040 


W4 


A34 


0.0025 


0.36 


0.73 


0.64 


0.27 


0.10 


0.006 


W4 


A 1234 


0.0027 


0.73 


0.98 


0.65 


0.10 


0.35 


0.021 



"^Difference maps from linear combinations of cliannels within a single DA, defined in equation (27). 
''Power in difference map in excess of white noise, \ ^Gi)t^^^^ — (C;)7gQ_^QgQ|. 
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Table 7. WMAP Boresight Pointing Vectors* 



DA/ bide 




rix 




Uy 




nz 


KIA 


0, 


.0399374 


0. 


.9244827 


-0. 


.3791264 


KalA 


-0, 


.0383635 


0, 


.9254372 


-0, 


.3769539 


QIA 


-0 


.0315719 


0, 


.9521927 


-0, 


.3038624 


Q2A 


0, 


.0319339 


0, 


.9522016 


-0, 


.3037965 


VIA 


-0, 


.0331733 


0, 


.9415643 


-0, 


.3351958 


V2A 


0, 


.0333767 


0, 


.9414947 


-0, 


,3353711 


WIA 


-0, 


.0091894 


0, 


.9394385 


-0, 


.3425944 


W2A 


-0, 


.0095070 


0, 


.9458644 


-0, 


.3244228 


W3A 


0, 


.0098004 


0, 


.9457678 


-0, 


.3246956 


W4A 





.0098081 


0, 


.9393480 


-0, 


.3428252 


KIB 


0, 


.0379408 


-0, 


.9239176 


-0, 


.3807057 


KalB 


-0, 


.0400217 


-0, 


.9246344 


-0, 


.3787473 


QIB 


-0. 


.0334030 


-0. 


.9517688 


-0, 


.3049925 


Q2B 





.0301434 


-0, 


.9519277 


-0, 


.3048361 


VIB 


-0, 


.0350363 


-0, 


.9409454 


-0, 


.3367405 


V2B 


0. 


.0314445 


-0, 


.9411385 


-0, 


.3365553 


WIB 


-0, 


.0114732 


-0, 


.9388325 


-0, 


.3441830 


W2B 


-0, 


.0115900 


-0, 


.9453501 


-0, 


.3258511 


W3B 


0, 


.0076818 


-0, 


.9454070 


-0, 


.3258014 


W4B 


0, 


.0075141 


-0, 


.9388923 


-0, 


.3441291 



''Beam line-of-sight unit vectors in spacecraft coordi- 
nates. Available in full precision in the released time- 
ordered data. 
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Table 8. Measured Gain and Baseline Susceptibilities'^ 





ag/9TFPA 


9b/9TFPA 




Radiometer 


(du mK-i) K-i 


mK K-i 




Kll 


-0.0021 


3.52 


0.1 


K12 


-0.0185 


5.05 


3.1 


Kail 


-0.0024 


-1.47 


0.2 


Kal2 


0.0077 


2.00 


-3.2 


Qll 


-0.0037 


3.79 


-1.1 


Q12 


-0.0016 


-3.52 


-1.6 


Q21 


0.0086 


-1.00 


-2.1 


Q22 


0.0058 


-0.57 


-4.6 


Vll 


0.0018 


57.4 


32.9 


V12 


-0.0045 


-6.23 


17.2 


V21 


0.0029 


6.10 


3.4 


V22 


-0.0002 


-9.43 


-3.8 


Wll 


0.0004 


-14.7 


-4.3 


W12 


0.00002 


-61.9 


-11.3 


W21 


0.0007 


-127. 


-7.4 


W22 


-0.0004 


-58.1 


0.4 


W31 


0.0003 


4.49 


5.3 


W32 


0.0021 


-20.2 


19.0 


W41 


0.0006 


41.4 


16.7 


W42 


-0.0011 


19.9 


8.0 



'^The thermal values are based on fits to a 10-day cooling period 
following a partial battery cell failure. See Figure 14 and §3.4.1. 



-57- 



Table 9. Limits on Spin-Synchronous Environmental Effects'^ 





Gciiii 


J_ lJ.l...i llLOjl. 


V UIUCL&C 


Band 


nK 


nK 


nK 


Kll 


-1.2 


22 


0.3 


K12 


-11.1 


32 


9.3 


Kail 


-4.2 


-9 


0.5 


Kal2 


5.8 


13 


-9.7 


Qll 


-1.3 


24 


-3.3 


Q12 


-1.2 


-23 


-4.7 


Q21 


31.0 


-6 


-6.3 


Q22 


35.2 


-4 


-13.8 


Vll 


0.7 


367 


98.8 


V12 


-6.0 


-40 


51.7 


V21 


5.9 


39 


10.3 


V22 


-1.7 


-60 


-11.3 


Wll 


4.0 


-94 


-13.0 


W12 


-0.1 


-396 


-33.9 


W21 


9.1 


-812 


-22.3 


W22 


-5.8 


-372 


1.1 


W31 


1.1 


29 


16.0 


W32 


1.1 


-129 


57.1 


W41 


21.1 


265 


50.1 


W42 


-30.1 


128 


24.0 


K 


6.2 


27 


4.8 


Ka 


0.8 


2 


4.6 


Q 


15.9 


2 


7.0 


V 


0.3 


77 


37.4 


w 


0.1 


17:-') 


9.9 



^la upper limits derived from measured gain and 
baseline susceptibilities in Table 8, combined with 
upper limits on temperature and voltage fluctua- 
tions at the spin pc;riod. Sign is preserved for each 
radiometer for roll-up by band. 
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Table 10. Upper Limits on Radiometer Cross Talk 



DA 


Electrical 


Radiometric 




dB 


dB 


Kl 


-37.7 


-26.8 


Kal 


-39.5 


-30.4 


Ql 


-41.6 


-32.3 


Q2 


-41.5 


-32.2 


VI 


-43.1 


-35.2 


V2 


-42.8 


-35.4 


Wl 


-48.8 


-48.3 


W2 


-47.1 


-43.5 


W3 


-38.6 


-42.6 


W4 


-46.1 


-47.5 



